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ABSTRACT 

The bimodality in observed present-day galaxy colours has long been a challenge for hi- 
erarchical galaxy formation models, as it requires some physical process to quench (and 
keep quenched) star formation in massive galaxies. Here we examine phenomenological 
models of quenching by post-processing the star formation histories of galaxies from 
cosmological hydrodynamic simulations that reproduce observations of star-forming 
galaxies reasonably well. We consider recipes for quenching based on major mergers, 
halo mass thresholds, gas temperature thresholds, and variants thereof. We compare 
the resulting simulated star formation histories to observed g ~ r colour-magnitude 
diagrams and red and blue luminosity functions from SDSS. The merger and halo 
mass quenching scenarios each yield a distinct red sequence and blue cloud of galaxies 
that are in broad agreement with data, albeit only under rather extreme assumptions. 
In detail, however, the simulated red sequence slope and amplitude in both scenarios 
are somewhat discrepant, perhaps traceable to low metallicities in simulated galax- 
ies. Merger quenching produces more massive blue galaxies, earlier quenching, and 
more frosting of young stars; comparing to relevant data tends to favor merger over 
halo mass quenching. Although physically-motivated quenching models can produce 
a red sequence, interesting generic discrepancies remain that indicate that additional 
physics is required to reproduce the star formation and enrichment histories of red 
and dead galaxies. 

Key words: galaxies:evolution - galaxies: luminosity function 



1 INTRODUCTION 



Galaxies in the local universe broadly fall into two main 
categories: blue star-forming disks typically found in lower- 
density environments, and red and "dead" ellipticals found 
in galaxy clusters. Some of the earliest extragalactic re- 
search revealed the two main morphological categories (e.g. 
iHubble 19261. and later work uncovered relations hips be- 
tween morphologies and lar g er scale environ ments (|Oemlerl 
1 1974 iDavis fc GeiieJ 1 19761 : iDressleJ Il980t ). More recent 
work employing data from large surveys has solidified a 
picture of galaxy bimodality, showing distinct divisions 
in galaxy colour (.Strateva et al., ,200 L : iBaldrv et al.. ,2004 : 
iBalogh et a"Lll2004l '). which serves as a proxy for star forma- 
tion rate per unit stellar mass. Surveys such as the Sloan 
Digital Sky Survey (SDSS) have characterized the statis- 
tics of red and blue galaxies to remarkable precision. Red 
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spheroids dominate the massive galaxy population (with 
stellar masses above ~ 10^'^'^ Mp,) , while less massive galax - 
ies tend to be star- forming disks l|Kauffmann et al]|2003d ). 
with luminosity functions showing that the relative fraction 
of red galaxies grows with luminosity. Surveys pushing to 
higher redshifts have robustly demonstrat ed that ga l axy bi- 
modality exi s ts ou t to redshift z = 1 (|Bell et al.l l2004bl : 
IWeiner et al] l2005l : IWillmer et all l2006t), and possibly to 
z = 2 and beyond (|Kriek et al.l l2008l : (Taylor et all l2009l : 
iBrammer et al.ll2009l ). The properties and evolution of mas- 
sive galaxies is one of the best- studied observatio nal areas 
of extragalactic astronomy (e.g. iFaber et al.ll2007l '). 



1.1 Why Is Quenching Needed? 

Despite the accumulating wealth of data, the physical origin 
of the bimodality in galaxy properties remains poorly under- 
stood. Red passive galaxies evolve from blue star-forming 
ones, in the sense that galaxies dominated by old stars to- 
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day must have built up their steUar mass through star for- 
mation at early epochs. Early theories of galaxy formation 
proposed that primordial clouds of gas collapsed under grav- 
ity to form rotati ng discs that harbor star formation (e.g. 
lEggen et al.lll96"3 ). In hierarchical mo dels, smaller galajcie s 
then merge to build up larger ones fe.g. lWhite fc Reesll978l ). 
Such mergers are thought to be capable of transforming spi- 
rals int o ellipticals through well- understo od dynamical pro - 
cesses l|Toomre fc ToomrdflQT^ : iMihos fc Hernauist|[l99^ . 
but it is less clear why this also results in a colour transfor- 
mation from blue and star-forming to red and dead. Recent 
efforts to measure the global mass build-up of the red galaxy 
population over cosmic time strongly suggest this type of 
transformation. The total stellar mass in red galaxies has 
grown by roughly a factor o f two since z~l, whereas blue 
galax ies have changed little dSell e t al] l2004bl : iFaber et al.l 
|2003). Furthermore, evidence s uggests the transition from 
blue to red is rapid, ~ 1 Gyr (jBell et all l2004bl : iBlantonI 
I2OO6I ). 

These data call for some mechanism(s) that "quench," 
or quickly shut off, star formation. Given the correlation 
between morphology and colour, such mechanisms (or re- 
lated ones) must also lead to morphological transformation 
of disk-dominated galaxies into bulge-dominated ones, ei- 
ther through disruption or disk fading and bulge growth. 
Requiring that the star formation quenching occurs rapidly 
likely excludes the possibility that galaxies simply consume 
all of their available fuel supply. Even without the rapid 
quenching requirement, analytic arguments suggest that the 
dense central regions of hot halos should cool rapidly onto 
galaxies in so-called cooling flows, which are not observed. 
Furthermore, simulations show that (barring new physics) 
filaments in the intergalactic medium can also continuously 
feed galaxies with fresh s upplies of cold gas (|Keres et al.l 
l2009bl : iBrooks et all 120091 '). This implies that feedback as- 
sociated with star formation is unlikely to quench galaxies, 
since red and dead galaxies have no star formation yet must 
remain quenched over the majority of cosmic time. Hence 
some additional physical process is required to quench star 
formation. 



1.2 Mechanisms that quench star formation 

There have been a number of mechanisms proposed to pro- 
vide quenching. Each mechanism either heats gas to the 
point where it cannot collapse to form stars ("preventa- 
tive" feedback) or expels gas that would otherwise form stars 
("ejective" feedback) (Kcrcs ct al. 2009a). In this work, we 
focus on quenching in massive central galaxies, so we do 
not address mechanisms that may quench star formation in 
small satellite galaxies su ch as ram pressure or tidal strip- 
ping and "strangulation" l|Gunn fc Gott|[l972l : iBalogh et al.l 
I2OOOI ). 

One well-studied mechanism is virial shock heating of 
intergalactic gas falling into a galactic dark matter halo. As 
gas falls in toward the central galaxy, its gravitational en- 
ergy converts to heat. In halos with masses above ~ 10^^ 
Mq, radiative cooling cannot keep up with gravitational 
heating, and a halo of hot virialized g as forms near the 
virial radius (iBirnbo im fc Dekell l2003l : iKeres et al.] l2005l : 
iDekel fc Biriib oim 200t^). This hot halo shocks infalling cold 
gas immediately to the viral temperature. Although it may 



succeed for ^2 Gyr, virial shock heating alone cannot quench 
all star formation, since the hot halo gas should cool via 
radiation a nd co ndense to form stars in the centres of mas- 
sive halos ((Birnboim et al.ii2007l ). Additional heating of the 
hot halo, perhaps by an active galactic nucleus (AGN; e.g. 
ICroton et al.ll2006|) or gravitatio n al heating due to c l umpy 
accretion (IBirnboim et all l2007l : iDekel fc BirnboimI 120081 : 
iDekel et al.ll2009bl ) could potentially prevent such cooling. 
This mechanism provides a physical explanation for the 
mass-dependence of the galaxy bim odality and the exponen- 
tial cutoff in lumi nosity functions l|Kauffmann et aLll2003cl : 
ISobral et al.ll201Gl ). 

As alternatives to virial shock heating, various types of 
galaxy interactions show promise for transforming galaxies. 
High-velocity encounters between galaxies may cause minor 
distortions and disk heating. Major galaxy mergers tend 
to transform disks to elli pticals (Toomrc & Toomrel 1 19721 : 
ISpringel et al.l l2005al : ICox et al.l l2006al ). and the resulting 
shocks and intense star formation winds may expel or heat 
cold star-for ming gas through shocks or feedback from su- 
pernovae (e.g. ICox et al.ll2006bl ) . Simulations show, however, 
that processes like supernova fe edback are not enoug h to 
sufficiently stifle star formation (|Springel et al.|[2005al ). As 
with virial shock heating, some additional energy input is 
required. 

A popular culprit to supply the necessary energy is 
the supermassive black hole thought to reside at the centre 
of most massive galax ies (cf. iKormendv fc Richstonelll995l : 
iRichstone et al.lll998l ). When the black hole accretes ma- 
terial, th e resulting AGN injects energy into its surround- 
ings (e.g. iMcNamara et ahlbOOd ) . supplying the nece ssary 
heating and/or expulsion of gas. ICroton et al.l l|2006l ) pro- 
pose that cooling and accretion of hot halo gas triggers low- 
luminosity AGN at the centres of massive elliptical galaxies, 
visible as radio lobes. This so-called "radio mode" feedback 
adds enough energy to the surrounding hot halo to prevent 

further cooling and star formation. 

In a me rger-driven que n ching scenario. iDi Matteo et ahl 

(2005) and iHopk ins et al. 1 l|200d ). among others, suggest 
that powerful AGN induced by galaxy mergers dump en- 
ergy into the remnants of the gaseous galactic disk, effec- 
tively superheating and expelling it to prevent further star 
formation. This "quasar mode" feedback operates only in 
the aftermath of major mergers of gas-rich galaxies. An ad- 
ditional benefit of this high accretion rate mode, for which 
the phenomenon was originally investigated, is that it ex- 
plains the existence of quasars and the relationships between 
black hole mass and bulge properties. 

Observational evidence provides some support for a con- 
nection between AGN and galaxy evolution. Correlations 
between black hole mass and bulge properties like mass and 
velocity dispersion helped motivate the concept o f AGN 
fuelin g and self-regulation via fe e dback |M agorri an et al.l 
19981: iFerrarese fc Merrit^ |200G| : iGebhardt et al.1 |200G| : 



Tremaine et all 120021 '). Furthermore, in both the local and 



distant universe, AGN have been ass ociated with a tran- 
sition in galaxies from blue to red (Sanchez et al.l |2004| : 
,Silverman e t al. 2008: Schawinski et al. 200 9) and from disk- 
dominated to bu l ge-dominated dKauffmann et al. I l2003al : 
iGrogin et al.ll2005l : iGabor et ahlbOOgf T 

A handful of alternative m echanisms have rec ently 
emerged that do not invoke AGN. IOekel et al.1 (|2009bl ') sug- 
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gest that smooth cosmic gas accretion provokes instabilities 
in a galactic disk by building up its mass, whereas clumpy 
accretion suppresses instabilities via gravitational stirring. 
Clumps transfer gravitational energy to the disk, effectively 
heating it and suppressing star form ation. Other m echa- 
nisms, such as cosmic ray heating (e.g. |Mathew&"2009|) and 
self-annhilation of dark matter particles (jTotani 2005), have 
not enjoyed as ample consideration in the recent literature. 



1.3 Hierarchical models of red and dead galaxies 

Several of these quenching mechanisms have been incor- 
porated into semi-analytic models (SAMs) of hierarchical 
galaxy formation. These models start with pure N-body sim- 
ulations of dark matter, then use the resulting dark mat- 
ter halo population and merger history to build baryonic 
galaxies within them using physically-motivated prescrip- 
tions. SAMs enable great speed and flexibility to explore 
parameter space and compare the resulting observables to 
data, but the large numbers of assumptions and free param- 
eters they employ make interpretation a challenge. 

Empirically-constrained SAMs have enjoyed success at 
reproducing obser vations of massive galaxies at low red- 
shift. The SAM of 'Croton et al.' ('200dl) based on the Mille- 
nium Simulation (Springcl et al. 2005b) established a red se- 
quence through mergers and disk instabilities, but required 
AGN feedback to prevent additional star-forming gas from 
cooling onto the most massive galaxies. This model success- 
fully reproduced luminosity functions derived from 2dFGRS 
ifMad gwick ct al. 200 2). Bower c t al. ( 200^ improved on 
the "Durham" SA M ((Cole et al.ll200d : iBenson et al.ll2002l : 
iBaugh et al.l l2005h , which forms bulges as the result of 
strong disk instabilities, by incorporating black hole feed- 
back in quasistatically cooling halos. Along with matching 
low-redshift luminosity functions and colour distributions, 
the improved Durham model matches observed luminos- 
ity functions, mass functions, and th e cosmic star forma- 
tion rate density at redshifts z < 2 JPozzetti et all l2003l : 
iDrorv et aLl l2003\ Cattanoo ct alj l|2006l ) improved on the 
GallCS SAM (Hatton ct al. 20o3) by imposing a critical 
halo mass above which accretion onto galaxies is shut down, 
and showed that this shutdown yields b etter matches to 
glaxjy color distributions and luminosities. ISomerville et al.] 
2008) incorporate a "unified" model of AGN feedback into 
their SAM, including both a low-accreti on "radio mode" an d 
high-accretion "bright mode" following ISiiacki et al1l|2007t ). 
This model matches physical properties of galaxies, such as 
the stellar mass function, along with black hole correlations 
with their hos t bulg es. Taking a more analytical approach, 
I Hopkins et al.l l|2008l ) adopt dark matter halo mass functions 
from simulations, populate halos with galaxies using halo 
occupation distribution models, and estimate merger rates 
among galaxies. By assuming that major mergers convert 
star-forming disk galaxies into spheroidal quiescent galax- 
ies, these models are found to match mass functions and the 
integrated stellar mass density of quenched galaxies. 

Overall, the general scenario t hat has enjoy e d the 
most success is as follows (cf. iHopkins et all l2008l : 
iKormendv et al]|2009l ): A major merger causes a transfor- 
mation from a spiral to bulge-dominated system, while con- 
currently feeding a central black hole that shines briefly as 
a quasar and emits a small fraction of its accreted mass- 



energy back into the galaxy. This causes a rapid truncation 
of star formation by driving a large fraction of star-forming 
gas from the galaxy, forming a red and dead elliptical. After 
the transformation, eventually a cooling flow forms, allow- 
ing low angular momentum gas to feed the black hole. This 
triggers a low-luminosity AGN which deposits energy into 
the surrounding hot gaseous halo. This shuts off the cool- 
ing flow, which in turn shuts off the AGN, and the cycle 
restarts. SAMs incorporating this scenario can explain ob- 
served correlations between galaxy morphology, colour, stel- 
lar age, bulge mass, and black hole mass. 

Recently, the growth of and feedback from black 
holes have been i ncorporated into hydrodynamic cosmolog- 
ical simulations (ISiiacki et ahllioOTl : iDi Matteo et al.ll2008l : 
iBooth fc Schavj 20091 ). At substantial computational ex- 
pense, these simulations combine with N-body dynamics the 
complex three-dimensional dynamics of gas inflow and out- 
flow responsible for establishing galaxy properties. Such sim- 
ulations employ parameters for subgrid physical processes 
associated with star formation, galactic winds, and feedback 
from black holes, but bypass many additional parameters for 
treating gas dynamics as are required in SAMs. 

Simulations incorporating AGN feedback have showed 
early promise by reproducing properti es of the black hole 
and AGN and galaxy populations ( Siiacki et al.l l2007l: 
Colbcrs & Di Matteo 2008: ' Croft et allboogl : iDeGraf et all 



2009 ; Booth fc Schave 2009.), along with the global cosmic 
star formation history (|Schave et al.ll2009l ). However, the im- 
plementation of AGN feedback is quite heuristic, and some 
assumptions such as spherical input of AGN feedback energy 
seem dubious at face value. Hence like with SAMs, it is pre- 
mature to regard these models as fully physical descriptions 
of black hole and galaxy co-evolution. 

In this paper, we take an approach that might be con- 
sidered as a hybrid between simulations and SAMs. We 
begin with galaxy star formation histories (SFHs) taken 
from cosmological hydrodynamic simulations, and then ap- 
ply quenching recipes by post-processing the SFHs to see 
which recipe(s) best match observations of massive galaxies. 
This retains the rapid computational flexibility of SAMs by 
allowing us to explore parameter space without re-running 
expensive simulations, while still employing star formation 
histories of galaxies that are described by fully hydrody- 
namic simulations up until quenching. We compare our sim- 
ulated quenched galaxies to observed colour-magnitude di- 
agrams and luminosity functions in the local universe, fo- 
cusing on the massive (and bright) galaxy population since 
these are well-observed, since the largest discrepancies in 
existing models are found there, and since it is where the 
bimodality is most pronounced. This approach avoids any 
explicit reference to the physical mechanisms of quenching, 
e.g. it does not explicitly account for growth and feedback 
from black holes, and instead focuses on asking the question: 
Which aspects of galaxy and halo evolution govern quench- 
ing in massive galaxies? 

Our paper is organized as follows. ^ describes our 
base cosmological hydrodynamic simulation employed in 
this work. ij3] details how we post-process the SFHs under 
various proposed quenching mechanisms. We then compare 
the resulting galaxy properties to data as described in J?) 
Our results, in ^ suggest that several quenching mecha- 
nisms can indeed qualitatively reproduce a red sequence. 
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but they all fail to exactly match observations in sometimes 
subtle ways. We devote much of our discussion in i]6.4l to un- 
derstanding the failures of these models, with an eye toward 
what new aspects may be required to improve agreement. 
Finally, we summarize in ij?) 



2 SIMULATIONS 
2.1 Input physics 

We run simulations using a modified version of the N-body + 
Sr aoothed Partic le Hydrodyn amics (SPH) code GADGET - 
2 (|Springelll2005l ) described in lOppenheimer fc Davj ||2008D . 
In essence, SPH is a Lagrangian implementation of hydro- 
dynamics where particles represent fluid elements, and each 
particle's density and temperature are determined by kernel 
averages over its nearest neighbors. GADGET-2 computes 
gravitational forces between particles using a tree-particle- 
mesh algorithm, and employs an entropy-conserving SPH 
formulation to compute pressures and densities of the gas. 
The code tracks three distinct types of particle — dark mat- 
ter, gas (or SPH), and star — including positions, velocities, 
and masses, as well as temperatures, densities, and metal- 
licities as appropriate. 

GADGET-2 incorporates a mechanism for gas parti- 
cles to spawn star particles stochastically. This mechanism, 
based on an analytic two-phase model of the interstellar 
medium with cold dense clouds emb edded in a diffuse ion- 
ized gas (|Springel fc Hernauistll2003l ). operates on scales not 
resolved by the simulations. The model assumes that each 
gas particle that is sufficiently dense to be Jeans unsta- 
ble contains a hot ambient medium, cold clouds, and stars 
which form within the cold clouds, and the phases exchange 
mass, energy, and metals vi a condensation and supernovae 
energy input (|McKee fc Ost rikcr 1977,). The resulting star 
formatio n rate (SFR) in the cold clouds is calculated as- 
summg a lSchmidtj (|l959l l law. A gas particle with a positive 
SFR will spawn a star particle with a mass, M* , with some 
probability determined by the SFR. Any newly spawned 
star particle inherits the position, velocity, and metallic- 
ity of its parent gas particle, whose mass is reduced by 
M». We use M« = 0.5Mgas, the initial mass of a gas par- 
ticle. On average, a gas particle can spawn up to two star 
particles, but this varies because gas particles can acquire 
addit ional gas mass from nearby stellar evolution (see be- 
low). ISpringel fc Hernguistl (|2003D have reduced this model 
to a single free parameter: the characteristic time-scale over 
which cold clouds convert into stars at the threshold den- 
sity, which is set to 2.1 Gyrs in orde r to match the observed 
Kennicutt relation l|Kennicuttlll998t ) . 

Galactic outflows driven by star formation are im- 
plemented in a si milarl y stochastic fashion. Our model, 
first introduced by Spxingel fc Hcrnquist (2003,), ejects star- 
forming particles from galaxies by kicking them with a ve- 
locity Vw The probability for ejection is set by the mass 
loading factor 77 (i.e. the mass loss rate relative to the 
star formation rate) times the star formation probabil- 
ity for that particle. The simulation used here employs 
our f avored momentum- driven wind scalings llMurrav et al.l 



(|2005h . In these scalings, Vw is proportional, and r] is in- 
versely proportional, to the galaxy circular velocity. We also 
decouple the winds hydrodynamically until they reach a den- 
sity 0.1 times the critical density for star formation (up 
to a maximum duration of 20 kpc/i;„,) in order to mock 
up chimneys through which outflowing gas can escape but 
that are poorly treated by the spherically-averaging SPH 
algorithm. Both outflows and star formation are to be re- 
garded as Monte Carlo prescriptions, in which individual 
events are not b y themselves meanin gful; resolution conver- 
gence tests (e.g. iFinlator et al]|2006l ) have shown that once 
64 star particles are formed within a galaxy, then its star 
formation history is fairly stable and well- resolved. We use 

this as our galaxy (stellar) mass reso lution limit. 

More recent improvements in lOppenheimer fc Pavel 

(|2008t ) include a sophisticated chemical enrichment model. 
Type la supernovae, and stellar mass loss. We track car- 
bon, oxygen, silicon, and iron individually, with yields taken 
from various works. Type la SNe rates are taken from 
IScannapieco fc BildstenI ((200^, with a prompt and delayed 
compenent, and these inject energy and metals (predomi- 
nantly iron) into surrounding ga s. We also track mass lo ss 
from AGB stars (using models of lBruzual fc Charlot|[2003l . ; 
BC03), injecting metal-enriched mass (but no energy) to the 
nearest three gas particles. Generally, these "delayed feed- 
back" mechanisms become more important at later epochs, 
which will be less of a concern for us in this work since the 
quenching of star formation in massive galaxies typically oc- 
curs relatively early on. 

The momentum-driven outflow scalings, together with 
hydrodynamic decoupling, have resulted in our simulations 
broadly matching a wide range of galaxy and intergalac- 
tic medium properties. These include the chemical enrich- 
ment of the 2^2 — 6 IGM (|Oppenheimer & Dave' '2OO1 
20081) the galaxy m ass-m etallicity relation (Dave ct al 



2OO7I : IFinlator fc Da vd 2008) , and high- 2 galaxy luminosity 



20051) as described i n lOppenheimer fc Davg (|2006l ) and 



Oppenheimer fc Davel (120081). an d mot ivated by observa- 



tions of local winds by iMartinI (|2005l ) and iRupke et al.l 



functions (|Dave et al.. ,2006. : iFinlator et al.ll2007h . As such, 
while our model is not an ab initio description of galactic 
outflows, it appears to modulate the SFHs of star-forming 
galaxies in broad accord with observations. Hence our sim- 
ulations provides a plausible starting point for studying the 
quenching of star formation across cosmic time. 



2.2 Simulation parameters 

The simulation we employ here begins at redshift 2 = 99 
with 512^ dark matter particles and 512^ gas particles in 
a cube of comoving side length Zbox = 96/i~^ Mpc with 
periodic boundary conditions, and is evolved to z = 0. 
The volume of this simulation allows us to sample large 
galaxies in high-density environments (with halo mass up 
to ~ 10^^'^Mq), while also providing the dynamic range 
to resolve galaxies well fainter than the knee of the lumi- 
nosity function. With a mass resolution of 1.2x10* Mq 
per gas particle, we can resolve galaxies with stellar masses 
^ 3.8 X lO^Mfi^. We em ploy a WMAP-5 concordance cos- 
mology (|Komatsu et all i2009) with Ho = lOO/i = 70 km 
s^^ Mpc~^, matter density Q.m = 0.28, baryon density 
Q,b = 0.046, a cosmological constant with Q.a = 0.72, 
root mean square mass fluctuation at separations of 8 Mpc 
(78 = 0.82, and a spectral index of n = 0.96. 
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2.3 Simulation outputs and analysis tools 

We configure GADGET-2 to output simulation "snapshots" 
at 108 redsiiifts, starting at 2: = 30 and ending at 2: = 0. 
The proper time between snapshots ranges from ~ 50 Myr 
at 2 ~ 6 to ~ 330 Myr at 2 ~ 0. The snapshots contain 
information for each dark matter, gas, and star particle in 
the simulation, including position, velocity, mass, metallicity 
(for gas and star particles), density and temperature (for gas 
particles), and time of formation (for star particles). We also 
store auxiliary information such as star formation rates and 
carbon, oxygen, silicon, and iron abundances. 

From this basic information we identify and ex- 
tract properties of dark matter halos and galaxies. 
We use SKID (Spline Kernel Interpolative DENMAX; 
|http:/ /www-hpcc. astro. washington.edu/tools/skid. html ) to 
identify gravitationally-bound groups of star-forming gas 
and star particles that comprise a galaxy. Halos are identi- 
fied using a spherical overdensity (SO) algorithm. Beginning 
at the most bound particles within SKID-identified galax- 
ies, halos are expanded radially until they reach the virial 
overdensity threshold. Smaller halos whose centres lie within 
the virial radius of larger ones are subsumed, while halos 
whose virial regions overlap but whose centres lie outside 
each others' virial regions divvy up the particles according 
which centre is closer as scaled by their virial radii. In the 
end, each halo has at least one galaxy by construction, and 
each particle belongs to not more than one halo. From these 
identified galaxies and halos, we calculate quantities such as 
masses, star formation rates, and metallicities by summing 
over all member particles. 

The star particles within a SKID galaxy are tagged with 
formation time and metallicity, and hence provide a series 
of single stellar populations (SSPs) from which one can con- 
struct the star formation and chemical enrichment history of 
each galaxy. This allows us to calculate galaxy spectra using 
the high- resolution grid in the stel lar population synthesis 
models of lBruzual fc Chariot! (|2003l ). We then measure pho- 
tometric magnitudes by convolving the galaxy spectra with 
broad-band filter curves. The code optionally computes ex- 
tinction due to dust in the galaxy; we explore the effects of 
different dust models in tj5.5l Since we mainly focus on red 
sequence galaxies where dust extinction is negligible (e.g. 
iLauer et al.ll2005l ). our dust prescription does not strongly 
impact our results. 



3 POST-PROCESSED QUENCHING 

For each simulated galaxy we compile a star formation his- 
tory to 2 = 0. Hence we can test models for quenching by 
modifying that SFH in post-processing, based on some set 
of prescribed rules. This approach allows us to assess to 
first order the feasibility of a wide range of quenching mech- 
anisms in a short amount of time, without the additional 
computational costs associated with building those mecha- 
nisms directly into the simulations. To describe our quench- 
ing models in detail, we first present a general explanation 
of the quenching process, and then give examples of the im- 
plementation for each quenching mechanism. 

Since we determine properties like galaxy stellar masses 
and broad-band luminosities from groups of simulated star 



particles, we can mimic the effects of quenching simply by 
removing star particles from consideration when extracting 
these bulk properties. Thus, at each snapshot, we examine 
the conditions under which each new star particle formed. 
If these conditions match the quenching condition, then we 
flag that star particle to indicate that it never should have 
formed. With this flag, we track the quenched star particle 
through all subsequent time steps. Finally, when we extract 
a galaxy's stellar mass and luminosity (at any particular 
output time step), we ignore any star particles which have 
been flagged as quenched. 

This method has the obvious disadvantage that any gas 
which forms into a quenched star particle will be locked 
up in that phantom star particle (i.e. a star particle that 
never should have formed), rather than remaining in its 
gaseous state to be tracked dynamically by the simulation. 
Furthermore, the additional feedback energy associated with 
quenching could have an impact on the surrounding gas that 
subsequently forms stars, which our post-processing tech- 
nique cannot account for. Physically, one can view this ma- 
terial as having been heated or expelled in such a way so that 
it is unable to return to a cool star-forming state within a 
Hubble time, and also does not interact signiflcantly with 
other infalling gas or galaxies. Clearly this cannot be true in 
detail, but this approach can still give useful first-order in- 
sights. Hence while this method has the advantage of compu- 
tational flexibility and speed, a fully self-consistent quench- 
ing mechanism must be incorporated into the simulations 
dynamically in order to properly assess all effects; this is 
under development. 

We study quenching mechanisms that have been re- 
cently considered in the literature. Namely, we specify con- 
ditions under which star formation should be quenched us- 
ing three different models: (i) quenching induced by galaxy 
mergers, (ii) quenching by virial shock heating in massive 
dark matter halos, and (iii) quenching due to the inability 
of gas shocked above a critical temperature and/or expelled 
in a wind to cool. For each of these mechanisms, we define 
a quenching condition based on an event or property that 
we can extract from the simulations. Then we ignore star 
particles whose formation satisfies the quenching condition, 
as described above. We refer to this as quenching the star 
particle. 

3.1 Quenching via galaxy mergers 

As mentioned above, various works have suggested that feed- 
back from intense star formation and AGN activity induced 
by galaxy mergers can heat and expel the cold gas from 
merger remnants, rapidly halting star formation. We apply 
this hypothesis to the simulations by quenching star forma- 
tion in galaxy merger remnants. To do this, we first iden- 
tify remnants of major mergers by tracing growth in galaxy 
stellar mass. For each galaxy in a given output time step 
tc, we identify its most massive progenitor galaxy in the 
most recent time step tp [tc > tp, and these times represent 
the age of the universe at the corresponding time step). If 
the galaxy and its progenitor have stellar masses M« and 
Mp, respectively, then a major merger occurred in the last 
timestep if M« = Mp{l + 1/r), where r is the critical mass 
ratio that separates major and minor mergers. We use the 
typical value r = 3 (jPasvra et al.ll2006l : I Woods et ai]|2006l : 



6 Gabor et al. 



iHopk ins et al]|2008l ) and a more extreme r = 4, and also 
require that the remnant be well-resolved in our simulation, 
with M* > 3.84 x IQ^Mq. Note that since galaxies also 
grow by star formation, this merging criterion can be con- 
sidered conservative in the sense that the actual merger ratio 
is probably higher (i.e. r is larger) than the assumed value. 

Motivated by recent work suggesting that rem- 
nants of galaxy mergers may re - form into star-formin, 
disks llSpringel fc Hernguistj l2005l: [Robertson et al.l 1200' 
iGovernato et al.l |2007| . l2009l : iHopkins et al.l l2009l 'l. we in- 
corporate an additional parameter, /gas- We assume that 
any merger remnant with a gas fraction Mgas/(Mgas -I- 
Afstars) > /gas rc-forms into a disk galaxy, and therefore is 
not quenched. If a galaxy meets the quenching conditions, 
then we set its most recent merger time, tmorgc = tp. We 



track tn 



for all galaxies through the evolution of the 



simulation. At any given time step, galaxies inherit fmcrgc 
from their most massive progenitors of the previous time 
step. Note that in some cases, merger remnants that are ab- 
sorbed by larger galaxies can effectively lose their remnant 
status; this is fairly rare except at the earliest epochs. 

After identifying merger remnant galaxies, we quench 
any star particle that forms in such a galaxy at a time 
iform > imerge Siucc the simulatious track the time of for- 
mation tform for each star particle and our tools match star 
particles to galaxies, we can easily identify those star par- 
ticles to be quenched. Our model assumes that, barring a 
subsequent merger with a larger galaxy, a merger remnant 
will never recover from whatever quenching process occurs. 
That is, we maximally halt all future star formation, ig- 
noring the possibility that cold gas later accretes onto the 
remnant to rejuvenate star formation. A more physically ac- 
curate model might account for the initial heating of gas in 
the galaxy during the merger, then permit coo ling as part of 
the s ubsequent evolution. Some authors (e.g. ICroton et al.l 
I2OO6I ) have proposed that heating due to an intermittent 
AGN prevents later cooling and star formation; we implic- 
itly assume that such a process operates to keeps merger 
remnants quenched. 



3.2 Halo Mass Quenching 

Some authors ( Birnboim et al.] l2007l : iDekel fc BirnboimI 
l2008l : lDekel et al 1 12009b')" have proposed that hot gas in mas- 
sive dark matter halos can halt the accretion of gas onto 
galaxies, cutting off the fuel for star formation. As gas falls 
toward the centre of the dark matter potential well, gravi- 
tational energy converts to thermal energy, heating the gas 
to temperatures near the virial temperature. In halos with 
masses Mu < 10^^ ''~^^Mq, rapid cooling prevent the virial 
shock from being supported, but in halos above that mass a 
stable ho t gas envelope can form (|Birnboim fc Dekelll2003l : 
iKeres et^ al. 2005, 2009b). The virial shock rapidly heats any 
further infalling gas, stifling accretion of cold gas onto galax- 
ies and thus quenching star formation. This mechanism ef- 
fectively quenches star formation in dark matter halos with 



Mh > Mc, where Mc - IO^^Mq. 

The X-ray emission from intracluster gas around mas- 
sive galaxies indicate that the gas is cooling at rates of tens 
or hundreds of Mq per year. However, young stars and 
reservoirs of cool gas where the cooling flow might be de- 
posited have not been observed (the cooling flow problem; 



see IPeterson fc FabianI |2006| . for review) , suggesting that 
an additional long-term heating mechanism acts near the 
halo centre. Low-luminosity AGN could in principle pro- 
vide this heating (Croton ct al. 2006. ). but the details of the 
heating process are poorly understood. Some works sug- 
gest that AGN jets generate pressure waves or magnetic 
fields within intracluster ga s to isotropize energy input (e.g. 
iRuszkowski et al.l l2004a b: Briigg en et~aLl bOOsI ). though it 
is not clear that such processes can operate in more typical- 

size d halos. 

iBirnboim et al.1 (|2007| ) suggest a scenario in which dark 
matter halos that reach Mc abruptly form a rapidly expand- 
ing shock, heating the supply of infalling gas and quenching 
star formation for ~ 2 Gyrs. Eventually the shock slows, the 
halo stabilizes, and ~ 10^^ Mq of gas quickly cools and ini- 
tiates a starburst, followed by a long-term quiescence. This 
can quench star formation in galaxies residing in groups, but 
galaxies in massive clusters require additional heating from 
AGN or clumpy accretion to prevent recurrent star forma- 
tion. 

In this paper, we explicitly avoid reference to any par- 
ticular physics of halo quenching, but are rather interested 
in generally testing whether quenching based on some halo 
mass threshold can match observations of passive galaxies. 
We therefore adopt a simple approach in which no stars can 
form in halos above Mc . To model this process in the simula- 
tions, we first identify all dark matter halos with Mh > Mc 
in a given output time step. Then we quench any star par- 
ticle which formed in one of those halos within the most 
recent time step. Again, since the simulations and analysis 
tools track the time of formation and corresponding dark 
matter halo of each star particle, we can straightforwardly 
identify star particles to be quenched. We treat Mc as a 
free parameter when comparing with observations, testing 
plausible values in the range Mc = 10^^-^~^^-^ Mq. 

This procedure implicitly quenches star formation in all 
galaxies in the dark matter halo, including both dominant 
(or central) and satellite galaxies. Therefore, when a small 
galaxy falls into a massive halo, we immediately truncate 
its star formation rather than allowing a slow decline lik e 
that expected for "strangulation" (e.g. ISimha et al.l 120091 ). 
Although the abruptness of this tr uncation does not r eflect 
realistic galaxy infall to clusters (|Balogh et al.ll2000l ). we 
are primarily interested in the population of massive cen- 
tral galaxies that dominate the bright end of the red se- 
quence today. We discuss variants of this mechanism, in- 
cluding quenching only central galaxies, in ij5.4l 

3.3 Quenching of hot mode and wind mode 

Galaxies in simulations show two main paths in density- 
temperature space f or acquiring gas t o fuel star forma- 
tion (|Binnevl Il977l : iBirnboim & Dcke] l2003l : [k ere s etai] 
' 200f^ : Hot mode, in which gas heats to near the virial tem- 
perature of the halo before cooling to form stars; and cold 
mode, in which gas never gets significantly above ^ 10^ K 
and radiates most of its gravitational pot ential en ergy in Ly- 
man alph a emission tFardal et al.. ,200ll : lYang et al. .20061). 
Recently, iQppenheimer et al.l (|2009bl ~ introduced recycled 
wind accretion, or "wind mode" , as a third accretion path. 
Wind mode constitutes gas that was previously ejected from 
a galaxy, and has been re-accreted to form stars, and in our 
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Figure 1. Distribution of Tmax for simulated star particles in 
galaxies with stellar masses in three mass bins: 10^"— lO'^'^'^ , lO'^^ - 
10^^'^, 1O^^-1O^^-^M0. All masses show two distinct peaks, with 
a low shoulder or tail to high temperatures. Vertical dashed lines 
show our two separate choices of the critical Tmax above which 
star particles will be quenched. 



simulations with plausible outflows models it dominates the 
cosmic accretion budget, particularly in massive galaxies, at 
-2 < 1- 

In either the hot mode or wind mode cases, it is 
conceivable that the gas should never re-accrete. This 
could be owing to some heat so urce that ke e ps ho t ha- 
los energized, as hypoth esized in iKeres et al.l (|2005l ) and 
iDekel fc BirnboimI l|2006h . or else it could be because there 
are numerical difhcultics in handling two-phase media and 
cold clumps in SPH (Agcrt z ct al. 200Z) by wh ich accretion 
is over-estimated. Indeed, iKeres et al] (|2009tJ ) found that 
the amount of hot mode accretion in entropy-conserving 
SPH is far less than in previous version of SPH, highlighting 
the nu merical uncertainties involved. lOppenheimer et al.l 
l|2009bh also speculates that wind mode should never re- 
accrete, again owing to numerical issues of cold clumps mov- 
ing through hot halo gas. It is conceivable that, since hot 
mode and wind mode are preferentially more important in 
larger systems, that simply removing such accretion modes 
could yield a bimodality with a clear red sequence. 

Our simulations track a parameter Tmax corresponding 
to the highest temperature achieved by that particle in its 
history before entering a galaxy. The simulations stop up- 
dating Tmax for gas particles ejected from galaxies via winds, 
even if they are heated above the previous maximum value. 
When a gas particle stochastically spawns a star particle, 
the star particle inherits its parent's Tmax at the time of 
spawning, retaining it for the remainder of the simulation. 
We thus track the maximum temperature achieved by the 
gas which formed any given star particle. To quench hot 
mode, we simply assume that any gas particle heated above 
some critical temperature, Tc, can never cool to form stars. 
For this mechanism, we simply identify and quench any star 
particles with Tmax > Tc. 

Figure [1] shows the distribution of Tmax for all star 
particles in galaxies in three different stellar mass bins: 
M, = 10^°, lO^S IO^^Mq. Each distribution shows two dis- 



tinct peaks (at log T ~ 4.4 and 4.9) and a tail to higher tem- 
peratures. The low-temperature peak is from gas that has 
never been substantially heated, and the high-temperature 
peak comes from warmer gas that then cools to form stars. 
The peaks are associated wit h hydrogen cooling and he - 
lium cooling, respectively (cf. [Sutherland fc Dopital Il993l ). 
Mild shocks during accretion likely heat the warm com- 
ponent, while the colder component s omehow avo i ds such 
shocks. Somewhat counter-intuitively, iKeres et al.l ()2009bh 
found that both low and high mass galaxies form predomi- 
nantly through cold mode, and that intermediate mass sys- 
tems (Mbaryon « 5 X 10^^ Mq) have the highest hot mode 
fraction. In the high mass case, it is because large galax- 
ies assemble from lower-mass galaxies that formed early on 
mainly via cold mode accretion. Based on this plot, we ex- 
plore critical values of log Tc = 4.7,5.4 (shown as vertical 
dashed lines), though it turns out that the results are not 
very sensitive to this choice. 

Since most of the high-temperature gas is heated via 
virial shocks, one might expect these results to mimic those 
of the virial shock heating prescription in i]3.2l In many ha- 
los, however, cold flows of gas from the IGM (particularly 
at early epochs) can penetrate the hot envelope to feed the 
central galaxy directly, s o the quenching is not complete 
(jKeres et al.H2005l . [2009bh . Our results in indeed sug- 
gest that cold flows are a significant factor in ongoing star 
formation in massive galaxies. 

Quenching wind mode is also straightforward, since we 
track which gas particles have been ejected in a wind, and 
this information is passed on to spawned star particles (see 
lOppenheimer et al.l l200S)b'. for implementation details) . Cer- 
tainly, it is the case that at least some material that par- 
ticipated in an outflow (particularly at early epochs when 
outflows were most prominent) is likely to have fallen back 
into a galaxy at late times. Here, we make the most extreme 
assumption that none of it ever falls back. 



4 OBSERVATIONAL CONSTRAINTS 

The ultimate goal of this study (and others like it) is 
to build a model for galaxy evolution that matches ob- 
served global distributions of galaxy properties. Accord- 
ingly, we compare the results of our simulations to well- 
studied quantities of the low-redshift (2 < 0.1) galaxy pop- 
ulation, primarily based on measured luminosities. We fo- 
cus on colour-magnitude diagrams (CMDs) and luminos- 
ity functions (LFs). Of course, a litany of other observa- 
tions could provide additional detailed constraints, such 
as luminosity - and c olour-dependent correla t ion functions 



(Zehavi et al. 2005; Weinmann et al .1 I200a 'Coop er et al 
2006; Phl ep_s ct al 20 6: Cooper et all 12007: Coil et al 



i2008l : iBrown et al.ll2008l : Iwilliams et al.ll2009l; ICooper et al 



'2009'), pr operties of post-sta,rburst galaxies dZabludoff et al" 



1996; Qu intero et"aLl l2004l: [siake et al.ll2004l : iBalogh et al . 
(20051 : I Yang et al.l l2008l : IWild et all l2009l), the total star 



form a tion history of the un i verse (e.g. Hopkins fc BeacomI 
I2OO6I : [Thompson et al.ll200ll . I2006I'). the star formation in- 
tensity distribution of galax;ies (|Thompsonll2002l ). and other 
observations of galaxy properties at higher redshifts. Given 
the challenges that even low-redshift CMDs and LFs present 
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to our models, we defer more detailed comparisons to future 
work. 

We use t he Va lue-Added Galaxy Catalog (VAGC; 
iBlanton et~aLl l2005b() of the Sloan Digital Sky Survey 
l| Adelman-McCarthv et al.ll2008l : iPadmanabhan et aljlioosf ) 
for comparison to our simulations. The DR6 version of 
the VAGC includes k-corrected absolute magnitudes in the 
SDSS ugriz + JHK bands for more than 2.6 million galaxies, 
including a special low-redshift sample of ~170,000 galaxies 
with distances of 10 Mpc to 150 Mpc (redshifts 
roughly 0.003 to 0.05). Because the low-redshift sample's 
volume is roughly comparable to (though larger than) the 
volume of our simulations, we use it (rather than the full 
VAGC) to provide observational constraints on our mod- 
els. We convert the reported absolute magnitudes (which 
use ft = 1) to our preferred cosmology (with h = 0.70) 
with Mh=.7 = Mh^i + 51og(0.7). Colours are then straight- 
forward, and we compute lum inosity functions using the 
1/Knax method (jSchmidtl [l96i ) and the Knax values pre- 
sented in the VAGC. 

Because it effectively traces stellar mass, we plot r-band 
absolute magnitudes in most of our plots, and we use g — r 
colours for CMDs. We compare observations with results 
from our simulation snapshots ai z = 0.025, chosen to fall in 
the middle of the observed low-redshift VAGC range, though 
the simulation results are insignificantly different at z = 
0. We refer to redshift 0.025 as z « 0, or "low-redshift," 
throughout the remainder of the paper. 

Stellar masses for gala xies in SDSS (incl uding those in 
VAGC) were determined bv lKauffmann et al.l ([2003b), using 
template fits to spectra including the effects of dust extinc- 
tion. We cross-correlate these publicly available data with 
the VAGC catalog to obtain stellar masses for our compar- 
ison galaxy sample. 



5 RESULTS 

To begin, we illustrate the current state of galaxy formation 
in cosmological hydrodynamics simulations without quench- 
ing in Figure [2l This shows CMDs (top) and LFs (bottom) 
for our simulations compared with observations from SDSS. 
We separate blue galaxies from red using a solid line drawn 
in the CMD, and we show a best-fit to the SDSS red se- 
quence as a dot-dashed line (see below). Our mass resolu- 
tion produces a diagonal envelope in the bottom right of 
the CMD panel, corresponding roughly to r ~ —20. This 
cutoff affects the blue galaxy LP fainter than r ~ —20.5, 
although the red galaxy LF is not impacted within our plot- 
ting r ange, brighter than r = —20. We follow iFinlator et al.l 
l|2006l ) to estimate LF un certainties for our simulation us- 
ing t he jackknife method (|Luptonlll993 : IZehavi et aLll2002l : 
[Wein berg et al. 2004). We sample the volume of our sim- 
ulation eight times, each time excluding one octant of the 
simulation and calculating the LF for the remaining 7/8 
of the total volume. The variance in the LF of the eight 
sub-samples provides an estimate of the uncertainties due 
to poisson noise and cosmic variance. Luminosity functions 
of the SD SS galaxies have typi cal uncertainties at the 0.1 
dex level (|Blanton et al.ll2005a^ . smaller than the symbols 
we use in the plot. 

This simulation fails to produce massive red and dead 
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Figure 2. Colour-magnitude diagram in g — r vs. r (top panel) 
and r-band luminosity functions (bottom panel) for our origi- 
nal simulation without quenching, compared to data from the 
SDSS VAGC. The solid line separates blue and red galaxies in 
the simulations, and the dot-dashed line shows a best-fit to the 
SDSS red sequence (see i|5.1l l. Greyscale in top panel shows SDSS 
VAGC galaxy density scaled logarithmically, and colour-coded 
points show simulated galaxies. In bottom panel, symbols repre- 
sent SDSS VAGC data, split into blue (blue triangles) and red 
(red diamonds) galaxies. Lines show simulated luminosity func- 
tions for blue (blue dotted) and red (red dashed) galaxies. Our 
mass resolution cutoff corresponds to a diagonal envelope in the 
CMD at about r = —20, and the cutoff affects our luminosity 
functions fainter than r ~ —20.5 (vertical solid line). Without 
quenching, our simulation produces too many bright blue galax- 
ies, and almost no bright red ones. 

galaxies as observed. Almost all galaxies occupying the 
red sequence regime are low-mass satellites th at have been 
quen ched primarily via strangulation (see e.g. ISimha et al.l 
l2009t ). with just a few massive systems that are a red ex- 
tension of the blue cloud, without evidence for a distinct 
bimodality. 

We reiterate that this simulation includes strong galac- 
tic outfiows from star-forming galaxies. Clearly, even such 
fairly energetic feedback is insufficient to quench star forma- 
tion in massive galaxies. Galaxies almost always have sup- 
plies of gas to fuel new star formation. Galaxies accrete cold 
gas from the IGM through filaments of the cosmic web (cf. 
iKeres et al.ll2009bl ). and hot gas in galactic halos may cool 
onto the central galaxies. In our simulations, much of the 
accreting gas today is recycled from earlier galactic winds 
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expelled from the galaxy lOppenheimer et al.ll2009bl V This 
illustrates that some other physical process must quench star 
formation in massive galaxies. 

When calculating galaxy magnitudes, we neglect the ef- 
fects of dust to highlight the intrinsic red sequence, i.e. the 
red sequence made up solely of galaxies with old stellar pop- 
ulations. Red galaxies without ongoing star formation con- 
tain little enough dust that we can neglect dust effects on 
colours for the intrinsic red sequence. Although dust tends 
to redden star-forming galaxies from the blue cloud, obser- 
vations suggest onl y ~ 10-20% of the red sequence come s 
from such galaxies (|Bell et al. I l2004al : iBrammer et IlI|2009D . 
In i]5.5l we will explore dust effects on the population of star- 
forming blue cloud galaxies. 

In the following three sections, we present CMDs and 
LFs for merger-based quenching, halo mass quenching based 
on a critical halo mass, and accretion mode quenching. We 
focus on the intrinsic red sequence successfully generated 
in the first two models, exploring some properties of the 
red galaxy populations formed. Because of its sensitivity to 
dust, we defer discussion of the blue cloud to >i5.5l 

5.1 Merger quenching 

Figure [3] shows the results of our merger quenching pre- 
scription with a mass ratio of 3:1 {r — 3) with /gas = 0.6, 1 
(top two panels) and r = 4 with /gas = 1 (third panel). 
CMDs in all cases exhibit a distinct red sequence with a 
shallow but nonzero slope, and a tail of bright red galax- 
ies, qualitatively in agreement with observations. However, 
the LP for the /gas = 0.6 case shows a significant excess 
of very bright red galaxies, which is only partly mitigated 
in the /gas ~ 1 case. Note that these are quite high gas 
fractions: Disk galaxies like the Milky Way today have gas 
fractions well below 0.6, and even at high-z one has to go 
to M t <~ 10^*^ Mq t o get typical gas fractions as high as 
50% l|Erb et al.ll2006l ). Our results prefer never re-forming 
a star-forming disk once a major merger happens. Since 
hydrodynamic simulations without AGN clearly show that 
a disk re-forms even with more modest gas fractions (e.g. 
[Robertson et al.l l2006l : iGovernato et al.] |2009| ) , this implies 
that for this quenching mechanism to work, something must 
heat or eject all the gas such that it not only stops form- 
ing stars at that time, but prevents the re-formation of a 
gaseous disk. 

In the third panel we make the model even more ex- 
treme by quenching star formation after all 1:4 mergers, re- 
gardless of gas content. This has only a minor impact on 
the resulting CMD (relative the 1:3 case). It decreases the 
number of very bright galaxies, which mitigates the discrep- 
ancy there, at the expense of worsening the agreement for 
moderately bright galaxies. 

The LF shows a related problem that all merger models 
fail to produce a sufficiently sharp knee at a characteristic 
magnitude M* as seen in the data. The model can reason- 
ably match the number of galaxies at the low-luminosity 
end, but it over-produces bright galaxies and underproduces 
galaxies at ~ L* . This is a fairly generic problem in this sce- 
nario, present even in our most extreme case, arising because 
mergers occur at a wide range of masses and only weakly 
pick out a characteristic scale around ~ L* . 

Merger quenching (along with all our quenching mod- 
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Figure 3. Merger quenching colour-magnitude diagrams (top 3 
panels) and luminosity functions (bottom) for intrinsic (dust-free) 
red sequence galaxies. As indicated in the lower left of the panels, 
the CMDs correspond to merger quenching with a gas fraction 
thresolds of /gas = 0.6 and 1, above which merger remnants re- 
form into star-forming disks. In the third panel we count 1:4 r = 
4) mergers as major, quenching star formation in the remnant. 
The dot-dashed line shows the best-fit to the SDSS VAGC red 
sequence. In the bottom panel, we show luminosity functions for 
red galaxies for each of the three quenching models (lines), along 
with the true red sequence from the SDSS VAGC (diamonds). 
For clarity, we show representative error bars for one simulation 
only. All models overproduce bright galaxies, and underproduce 
fainter ones. 

els) exhibits what we call the "blueness problem": the red 
sequence is too blue by ~0.1 magnitudes in g — r colour. 
This corresponds to ~ 10% error in the flux ratio between 
the two bands. We explore several possible explanations for 
this problem in ^6.21 which likely has an origin in the over- 
all calibration of galaxy metallicities. Since it occurs in all 
quenching models, we cannot use this blueness problem to 
constrain quenching mechanisms. 

To more quantitatively compare the slope, scatter, and 
average colour, we perform linear fits to the simulated and 
real red sequences. We first separate the red sequence from 
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Figure 4. Slope and scatter of the red sequence, for SDSS data 
(circles) and two simulations: one with merger quenching with 
/gas = 1 and r = 3 (triangles), and another with halo mass 
quenching with Mc = 10^^ Mq (squares). Our models produce 
red sequences that have shallow slopes, g — r colours that are too 
blue, and smaller scatter compared to the true red sequence. 



the blue cloud using a constant g — r colour, 0.65 for real 
galaxies and 0.55 for simulated ones. Various reasonable 
choices of this separator, including those with nonzero slope 
chosen simply by eye, yield similar results. After dividing 
the absolute magnitudes into bins of width 0.5, we calculate 
the mean, median, and standard deviation of <; — r colour 
within each bin for galaxies above our separator. We fit the 
median points to a straight line, as shown in the top panel 
of Figure[4]for SDSS galaxies and our 3:1, /gas = 1 merger 
model. Along with highlighting the blueness problem, the 
figure shows that the simulated red sequence slope is some- 
what too shallow, at least in this quenching model. Note (in 
Figure [3| that the model with lower /gas produces a steeper 
slope, and a smaller blueness problem, albeit failing more 
spectacularly in the LF comparison. 

In the lower panel of Figure 21 we plot the scatter as 
a function of absolute magnitude. The observed scatter is 
~ 0.05 ma gnitudes, b u t part of that is observational un- 
certainties. I Cool et ahl (|2006l ) estimated that the intrinsic 
scatter in g — r is 0.035 mag, at least for the most lumi- 
nous (> 2.2L*) galaxies. This is still larger than the scat- 
ter in our simulations, which is good because our quench- 
ing prescription makes the maximal assumption of zero star 
formation since the time of quenching. Early-type galaxies 
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Figure 5. A closer look at the CMD red sequence resulting from 
3:1 mergers with /gas = 1. Galaxies are colour-coded by r-band 
luminosity- weighted mean stellar ages (top panel); tyoung, the 
age of the youngest (most recently formed) star particle (middle 
panel); and r-band luminosity- weighted mean stellar metallicities 
(bottom panel). For reference, solar metallicity is Zq ~ 0.012. 
Although the brightest red sequence galaxies generally have old 



stellar populations, they have traces of young stars {ty 



< 1 



Gyr). A distinct gradient in metallicity creates most of the spread 
in g — r colour. 



toda y do appear to have a small "fro sting" of star forma- 
tion (|Trager et al.ll2000l : lYi et al.ll2005l '). which would tend to 
increase the scatter. However, the simulated scatter is not 
significantly smaller than observed, so there is scant little 
room for additional scatter from frosting. 

Figure [S] explores the red sequence produced in the 3:1, 
/gas = 1 merger quenching model in greater detail. We se- 
lect this as our preferred merger quenching model since it 
yields the least deviations from the observed red galaxy LF. 
The top panel shows the r-band luminosity-weighted mean 
stellar age; the middle panel shows the simulated galax- 
ies in colour-coded b ins of iyoungi the age of the youngest 
star particle in the galaxy; and the bottom panel shows the 
r-band luminosity-weighted mean stellar metallicity. The 
metallicities shown are absolute mass fractions of elements 
heav ier than helium, where solar metallicity is Zq «0.012 
( Asp lund et"aLl |2005| ) . 

At the brightest end, the galaxies have old stellar 
populations on average. Less massive galaxies show a 
wider spread in age, and perhaps counter-intuitively, the 
intermediate-age galaxies are redder than the oldest galax- 
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ies. There are some small (i.e. ~ L*) galaxies that are fairly 
young, having just been quenched onto the red sequence. 
The middle panel shows that even the brightest galaxies 
generally have some quite young stars (< 1 Gyr old). This 
is because they live in the largest halos that are still assem- 
bling at the present day, and hence they recently subsumed 
galaxies that were recently forming stars (i.e. unquenched 
galaxies); we discuss this further in ^ In metallicity, the 
red sequence shows a gradient where the reddest galaxies 
are the most metal rich, and this occurs at all luminosities. 
This reflects the well-known fact from population synthesis 
models that metallicity is the primary determinant of colour 
in old stellar populations. Here we can see why intermedi- 
ate age systems are redder than the oldest systems: They 
generally have higher metallicities, having formed stars up 
until a later cosmic epoch. 



5.2 Halo mass quenching 

Figure [6] shows CMDs and LFs for halo mass quenching, 
analagous to Figure [3] for mergers. We show results for 
three values of the critical cutoff mass, 
and 10^^-^A'/q. As with mergers we obtain a distinct red 
sequence in qualitative agreement with data. The slope is 
slightly shallower than in the merger case, essentially zero, as 
shown in Figure [4] and the scatter around the red sequence 
drops to nearly zero for the brightest galaxies, in contrast to 
merger quenching where the scatter is independent of lumi- 
nosity. The blueness problem is present at roughly the same 
level as in merger quenching. 

In halo mass quenching, most red galaxies tend to 
clump in a relatively small range of absolute magnitude. 
For logMc=12.5 this clumping occurs around r = —22 (and 
stellar mass ^10^^ Mq), and scales with halo mass for the 
other Mc cases. Variations in star formation and merging 
histories smear out this clumping somewhat, but this gen- 
eral feature remains; no such feature is seen in the data. 
The bright end of the red sequence grows large through dry 
mergers after the quenching process, and faint galaxies are 
mostly satellites. The dark matter halos of central galaxies 
tend to achieve the critical mass when the galaxies have a 
particular stellar mass, and then move from the blue cloud 
to the red sequence. Thus, at all epochs, galaxies tend to 
move onto the red sequence at an effective critical stellar 
mass or absolute magnitude that corresponds to the critical 
halo mass. Variations in the time of quenching and metallic- 
ities then cause significant vertical scatter in g — r colour for 
these clump galaxies. At a fixed stellar mass, galaxies early 
in the universe (high redshif t) have lower metallicities than 
those later on (low redshift) (|Erb et al.ll2"00^ :[ Maiolino et al.l 
|2008| ) , so galaxies that move onto the red sequence later have 
higher metallicities. 

The LF is quite sensitive to M^: Mc = 10^^ Mq fits the 
luminosity functions best, with the other cases far underpre- 
dicting or overpredicting the number of bright galaxies. The 
CMD clump is evident in the LF as well. Even in the gener- 
ally best-fitting Mc — 10^^ Mq case there is a tail of bright 
galaxies that leads to an excess in the luminosity function 
at r < —23, just as in merger quenching. Overall, the lumi- 
nosity function produces a slightly steeper drop-off than the 
merger quenching models, although a characteristic knee is 
still less pronounced than in the data. 
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Figure 6. Halo mass quenching CMDs and LFs, analogous to 
Figure [3] We quench star formation in halos above the critical 



mass, Mc, for three values Mc 



10^ 



,10^ 



and IW^-^M, 



&■ 



Galaxies tend to enter the red sequence in a clump defined by a 
characteristic r-band absolute magnitude (or stellar mass) asso- 
ciated with Mc (e.g. r ~ -22 for Mc = 10^^-^). The luminos- 
ity functions (bottom panel) show an unobserved uptick at the 
brightest bins for all three cases. The model with Mc = IO^^Mq 
yields the best match to the SDSS luminosity function. 



Figure [7] shows the red sequence in more detail for the 
best-fitting log Mc = 12 case, with galaxies colour-coded by 
mean stellar age, iyoung, and mean stellar metallicity (from 
top to bottom). Like in merger quenching (cf. Figure [SJ, we 
find that the most massive galaxies have old stellar popula- 
tions. Unlike in merger quenching, these old massive galax- 
ies do not contain any young stars. Each of these galaxies, 
lying at the centre of a cluster, tends to merge only with 
smaller galaxies that have already been quenched a while 
ago because they live in the same dark matter halo above 
the critical halo mass. Thus, when the central galaxies grow 
through accretion of satellites, they do not obtain any stars 
younger than 1 Gyr. 

The metallicity panel shows that the highest metallic- 
ity systems are generally the smallest (and reddest) systems. 
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Figure 7. A closer look at the CMD red sequence resulting from 
halo mass quenching with a critical mass 

of lO^^A^Q, Colour cod- 
ing is the same as in figure [5] Unlike in our merger quenching 
model, the brightest galaxies include no young stars. 



The large, old systems typically have fairly low metallicities. 
Qualitatively, this explains why the red sequence slope is 
basically zero; the age and metallicity gradients cancel each 
other. The tight correlation between luminosity and metal- 
licity directly translates into a scatter in the red sequence 
that goes from nearly zero at the brightest end to fairly large 
at the faint end, as seen in Figure 21 



5.3 Accretion mode quenching 

Figure [8] shows results of Tmax (i.e. hot mode) and wind 
mode quenching. The top two panels show CMDs for two 
different threshold temperatures, 10* '^ K and lO^ * K. These 
are chosen as the bimodality separation in Tmax in our sim- 
ul ations and the (n o wind, no metal cooling) simulations 
ofH eres et al] l|2009bi" ) . respectively. The third panel shows 
the extreme case of quenching both hot mode with a lO'*'^ K 
threshold, and quenching all wind mode as well. 

Quenching hot mode alone does not produce a red se- 
quence; it is not bimodal, and the number densities of red 
galaxies are far too small at all but the largest masses. This 
is because cold mode, either pristine or in the form of re- 
cycled winds, continues to provide significant accretion at 
late times. Hence the idea that simply keeping gas hot in a 
hot halo reproduces a red sequence does not appear to be 
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Figure 8. Tmax CMDs and LPs, analagous to Figures |3] and 
[6] We quench star formation from gas particles that achieved a 
maximum temperature above a critical value, Tc, for Tc = 10*'^ 
and 10^ * K. These first models fail to produce a substantial red 
sequence because at all masses, galaxies accrete fresh cold gas 
from the intergalactic medium and form it into stars. Although 
this model may suppress global star formation, it fails to suppress 
all the star formation in a given galaxy. In the third panel we show 
a model where, along with lO''-'^ threshold, we also 

quench star formation from any gas particles that were recycled 
from a galactic wind. Although this model shows a reasonable 
red sequence in the CMD, it strongly suppresses the formation of 
blue cloud galaxies (not shown). 



viable, since it takes only a fairly small amount of ongoing 
star formation to make a galaxy blue. 

Quenching both hot and wind mode produces some- 
thing that looks like a red sequence, with an amplitude, 
slope, and scatter that is (coincidentally) comparable to the 
merger quenching case. However, the red galaxy LF contin- 
ues to have the wrong shape; it is roughly just a constant 
factor higher than i n the hot mode-only quenching case. 
lOppenheimer et all ()2009bl l anticipate this failure of wind- 
mode quenching, finding that allowing wind mode accretion 
leads to the best-matching stellar mass functions below the 
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turnover mass M *. iKeres et all (|2009lJ) highlighted a nu- 
merical problem in GADGET-2 in which cold clumps form 
owing to thermal instabilities that then rain down ballisti- 
cally onto the central galaxy; the trouble is that the clumps 
always occur near the resolution limit, even as the resolu- 
tion varies substantially. H ence this cold drizzle may be a 
numerical artifact, although iKeres fc Hernguistl (|2009l ) used 
very high-resolution simulations to show that at least some 
of it is likely to be real. 

Accretion mode quenching appears to be the least 
promising of our various quenching mechanisms, so we do 
not elaborate on these results any further. In the discussion 
we examine some implications of its failure ( i]6.4p . 

5.4 Variants on quenching models 

After considering our basic models with few free parameters, 
we attempted several variations on our quenching models to 
try to improve agreement with observations. In the end there 
were no obvious successes, though the results illustrate some 
interesting trends. 

As discussed above, halo mass quenching leads to char- 
acteristic clumping in the red sequence not observed in na- 
ture. This clumping arises from the sharp critical dark mat- 
ter halo mass cutoff. While observational uncertainties could 
smear out the clumping, we investigate two more physically- 
motivated reasons why the critical mass cutoff might not be 
so sharp. 

Physically, the halo mass threshold is representative 
of the mass scale above which a stable hot halo forms. 
Preventative feedback mechanisms can operate more ef- 
fectively in the presenc e of a hot halo (|Keres et al.l l2005l : 
iDekel fc Birnboimll2006l ). However, there is non-trivial scat- 
ter between halo mass and hot gas fraction (e.g. see 
iKeres et al.] l2009bh . Hence we tried a model variant in 
which we quench galaxies in halos with a hot gas fraction 
/hot = Mhot/(A^hot + Mco\d) above a critical fraction. We at- 
tempted models with critical fractions 0.2, 0.4, 0.6, and 0.8, 
which correspond roughly to halo masses lO^^'^-lO'^^ '^ Mq 
but with a scatter of ~0.2-0.3 dex. The results are shown in 
Figure [9] (with the CMD in the top panel), for a critical hot 
gcis fraction of 0.6. This variant smooths out the characteris- 
tic absolute magnitude or stellar mass at which blue galax- 
ies move to the red sequence, yielding a less clumped red 
sequence. However, the slope and scatter (and trends with 
luminosity) remain essentially unchanged. Also, the LF still 
shows the excess of very bright galaxies as noted in the halo 
mass quenching case. Hence besides the aesthetic appeal of 
removing the clump in the CMD, it does not fare any better 
(or worse) than halo mass quenching. 

As another variant, we examined a scenario in which we 
quench only in central galaxies of halos, not satellite galax- 
ies. This is shown in Figure [9] (second panel) for a critical 
mass of IO^'^Mq. The CMD in this model shows greater 
scatter in the red sequence and a much flatter luminosity 
function slope than our preferred model where satellites are 
quenched as well. Galaxies are able to build up more mass 
before ending up on the red sequence, and massive central 
galaxies end up merging with more massive satellites than in 
the preferred model. This failure is disappointing, as most 
scenarios of halo mass quenching only quench the central 
galaxy, and observations favor a more gradual quenching of 
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Figure 9. Red Galaxy CMDs and LPs for variants of our quench- 
ing models. The top panel shows the red sequence created by 
quenching star formation in haloes with a hot gas fraction above 
0.6, the second panel shows halo mass quenching where only cen- 
tral galaxies undergo quenching, and the third panel shows a 
model where both a merger and a halo mass above IO^^Mq are 
required to quench star formation. Although these models show 
promise at matching the LF at r > —22, they share the prob- 
lems of our simpler models at the bright end; an excess of bright 
galaxies. 

satellites fe.g. IWeinmann et al.l [20091 ). Indeed, SAMs that 
quench satellite galaxies in halos have diff iculties reproduc- 
ing the colour distribut ions of satellites (|Weinmann et all 
l2006l : lBaldrv et aLll2006l ). which is a drawback that is shared 
by our original halo mass quenching scenario. We do note 
that the red sequence slope seems to be in better agreement 
with data, as larger galaxies acquire larger (unquenched) 
satellites that are more metal-rich. 

We also investigated a hybrid model that requires both 
a major merger and a critical halo mass for quenching. That 
is, we quench only if a major merger occurs inside a halo 
with mass > lO^^M© (here we assume /gas ~ 1). Physically, 
this model mimics the action of a hot halo to prevent post- 
merger gas accretion from the IGM, perhaps by continued 
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Figure 10. Red galaxy CMD and LF for variable mass halo mass 
quenching. Instead of a constant critical halo mass above which 
shocks quench star formation, the crit ical mass varies wi t h red - 
shift. The top panel uses a model from lPekel fc BirnboinJ l|2008l ) 
where the critical mass required to sustain a hot halo increases 
with cosmic time. This model yields a dearth of galaxies fainter 
than r = —22. In the other two models, cold flows are assumed 
to penetrate hot haloes at high redshift due to the higher relative 
densities of cosmic filaments. This effectively causes the critical 
mass above which quenching occurs to decrease with cosmic time. 
We show two differe nt versions of the cr i tical mass as a function 
of red shift, based on iDekel fc BirnboimI l|200d) and lOcvirk" et al.l 
||2008| ). which differ only by the metallicity assumed for cosmic 
gas. These models exhibit a nearly-constant slope in the LF. 



energy injection from an AGN. In our original quenching 
model, such post-merger accretion is always prevented, but 
here we only prevent accretion in larger halos that can sup- 
port a virial shock. This model is perhaps closest to what 
is envisioned in current ideas for quasar plus radio mode 
quenching, in which a merger initially transforms the galaxy 
into an elliptical and quenches star formation, but only when 
the resulting halo is sufficiently large to form a virial shock 
does the resulting halo gas stay hot e.g. via low-level AGN 
activity. 



Figure (9] third panel, shows the results of this model. 
The results are fairly similar to halo mass quenching, show- 
ing that it is this aspect that is the limiting factor for quench- 
ing; mergers are frequent enough that it does not add a strin- 
gent criterion. The red sequence shows slightly more scatter 
than in halo mass quenching, and there is no evidence for 
clustering at a particular r-band magnitude; the slope is still 
incorrect, and (as always) the blueness problem persists. The 
LF is also similar to the halo mass quenching case, with an 
excess of very bright galaxies and a possible dearth of ~ L* 
systems. This model may, in fact, be the best-fitting model 
as well as the most physically-motivated, but the results are 
qualitatively very similar to those of our simpler models. 

Finally, we consider three variants of halo mass quench- 
ing where the critica l halo mass varies with r edshift. The 
first is motivated bv iDekel fc BirnboimI (|2008l ). who argue 
that the critical halo mass for quenching owing to gravita- 
tional clump heating varies slightly with redshift due to the 
evolving cosmic density of gas (see their Figure 3). In this 
scenario, the quenching mass varies from Mc — W^^'^Mq 
at 2 = 3 to Mc ~ IO^^ ^Mq at 2 = 0. For 2 > 3 we use the 
lower value, = IO^'^Mq. Figure IIOI top panel, shows 
the results. The shape of the LF is affected most strongly, as 
there is now a significant dearth of red galaxies fainter than 
r ~ —22. The high Mc at late times makes the fainter end of 
the red sequence similar to the fixed Mc = lO^^'^M© case, 
i.e., it does not populate the red sequence enough. We do 
note that it produces a larger red sequence slope that is in 
better agreement with data, although the scatter is larger. 
Overall, this model does not fare as well as our favored halo 
mass quenching model with a constant Mc — 10^'^ Mq, al- 
though it does illustrate that varying the quenching mass 
with epoch can produce notable changes in the red sequence 
(at the aesthetic expense of introducing more parameters). 



In another paper, iDekel et al.1 (12009^ 1 argue that Mc 
evolves with redshift in the opposite sense, namely that 
Mc ~ lO^^ '^AfQ out to z = 1.5, but at higher redshifts it 
increases rapidly owing to the ability of cold streams to pen- 
etrate t hrough hot halos at early epochs (see their Figure 5, 
and also lDekel fc BirnboimluOOfil ) . From an eyeball estimate 
of their conjectured hot halo mass limit in the presence of 
cold streams, we take Mc — IO^^ '^Mq for z < 1.5, and 
log Mc ~ 1.2z + 9 . 9 for z > 1.5. We also show a variation 
from lOcvirk et al.l HqoI) where gas in cosmic filaments in a 
hydr odynamic simulation is fou nd to have a lower metallicity 
than lDekel fc BirnboimI ^2006) assume. In this formulation, 
Mc ~ 10"-^ for z <3, and log Mc ~ lAz+7.3 for z > 3. The 
lower panels of Figure [10] show the results in these cases: A 
red sequence is produced, but the CMDs have a very small 
colour scatter, and more dramatically the LFs are power 
laws! The LFs of the two parametrizations differ primarily 
in their normalization. These models exacerbate the prob- 
lem of producing a knee in the LF over constant- Mc models. 

None of these quenching model variants stands out as 
obviously superior to the others. Since the simple models 
perform as well as the variants, and because they involve 
simpler prescriptions, we use merger quenching with fgas = 
1 and r — 3 and halo mass quenching with Mc = 10^^ Mq 
as our preferred models. For the remainder of the paper, we 
focus on these two models. 
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Figure 11. Dust comparison of CMDs and LFs for merger quenching with /gaa = 1- Our dust model (right panels) shifts bright blue 
galaxies up and to the right in the CMDs, introducing contaminants to the bright end of the red sequence. Arrows illustrate the change 
in position (from tail to head) due to dust for three galaxies, with star-formation rates fs3, 9, and 15 Mq y^~^ (right to left). The 
luminosity functions reflect the increase in bright red galaxies due to contamination, as well as the suppression of bright blue ones. 



5.5 The Blue Cloud 

So far we have focused on the red sequence. Our merger and 
halo mass quenching models also produce a star-forming 
blue cloud of galaxies that is separated from the red se- 
quence, as shown in the left panels of Figures [ll] and 1121 
Here we examine the properties of the blue cloud in these 
scenarios. 

Computing the observed luminosities of blue cloud 
galaxies is complicated by uncertainties in dust extinction. 
Our simulation tracks metallicity, which is correlated with 
dust extinction. There are also empirical correlations be- 
tween UV or blue-band luminosity and extinction. Since in 
general the bright blue galaxies are quite metal-rich, the 
exact form and nature of the assumed extinction makes 
a significant difference in their resulting brightness and 
colours. We only apply an extinction correction to blue 
gala xies, leaving qu e nched galaxies unaffected. Also, we em- 
ploy (CajzettFeFal] (|2000l ) reddening law where the optical 
depth due to dust varies as wavelength A""'^. 

We first tried employing a correlation between metallic- 
it y and extinction from SPSS, including scatter, as described 
in lFinlator et al.l l|2006l ). However, we found that the result- 
ing blue cloud had a very large scatter, and a substantial 
number of previously blue galaxies ended up redder than 
the red sequence. It is possible that second parameter cor- 
relations exist in metallicity versus extinction that are not 
accounted for, but we did not explore this further. 

Our preferred method, i.e. one that produced a blue 
cloud similar to observed, uses an empirical correlation 



between UV lum i nosity and dust extinc tion given by 
'Wa ng fc HeckmanI l|l996l ). also described in iFinlator et al.l 
(i2006r ). This prescription moves bright blue galaxies more 
than dimmer ones, so that the bluest galaxies are not nec- 
essarily the brightest ones. The resulting overall shape of 
the blue sequence matches the observed SDSS blue sequence 
reasonably well. 

Figure [TT] compares CMDs and LFs for our merger 
quenching model with and without dust. Without dust, the 
brightest blue galaxies in the universe would have r-band lu- 
minosities brighter than the brightest red galaxies, and they 
would be ^ 0.3 magnitudes bluer than the red sequence in 
g — r. With dust applied, the brightest star-forming galax- 
ies are scattered into the red sequence, making up ~6% of 
red sequence galaxies; this is reasonably consistent with ob- 
servational estimates (cf. iBell et aL l l2004al :lB rammer et al.l 
2009). Furthermore, the dust obscuration ensures that mas- 
sive blue galaxies are not brighter than the red sequence. 
The luminosity functions refiect this difference markedly. In 
the no-dust case, blue galaxies dominate over the bright end 
of the red sequence, but the presence of dust shifts the bright 
blue objects to lower luminosities. 

Figure [12] compares the dust and no-dust cases for halo 
mass quenching. Even the no-dust case lacks bright blue 
galaxies like those seen in the merger model and in ob- 
servations, and once dust is included the problem becomes 
significantly more severe. As discussed in >i5.2l once galax- 
ies attain a stellar mass corresponding to the critical halo 
mass, they move out of the blue cloud and onto the red 
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Figure 12. Dust comparison of CMDs and LFs for lialo mass quencliing witti Mc = IO^^Mq. Arrows stiow the cliange in position 
for three galaxies, with star-formation rates ~2, 4, and 7.5 Mq yr~^ (right to left). Even without dust (left panels), this model does 
not produce exceptionally bright blue galaxies: once blue galaxies attain a critical stellar mass corresponding to Mc, they move to the 
red sequence. Our dust model (right panels) suppresses the bright blue galaxies, leading to a steep cutoff in the blue galaxy luminosity 
function. 



sequence. Application of the dust prescription pushes blue 
cloud galaxies up and right in the figure, creating a diagonal 
envelope. Since dust tends to move the brightest galaxies the 
most, the dearth of blue galaxies intrinsically brighter than 
r ~ —22 leads to fewer intrinsically blue interlopers on the 
red sequence (< 2%). 

As an aside , we note an interesting study by 
iMaller et al] (|2009l ) of galaxy orientation in relation to dust 
obscuration, and its impact on derived galaxy properties. 
Notably, they find that structural parameters like axis ratio 
and Sersic index significantly impact measured colours and 
magnitudes, and they find an average of 0.2-0.3 magnitudes 
of extinction in the SDSS g and r bands. This study pre- 
dicts considerably more bright blue galaxies than are gen- 
erally inferred from studies that do not take into account 
orientation-dependent effects. We do not account for such 
effects in our work, as we use a simple dust screen model to 
account for obscuration. 

While such comparisons are illustrative, the sensitiv- 
ity to the exact extinction prescription used makes robust 
interpretation of the discrepancies difficult. As an alterna- 
tive approach, we can perform the comparison versus stel- 
lar mass functions rather than luminosity functions. The 
stellar masses derived by iKauffmann et al] (l2003a ) implic- 
itly account for dust extinction on a galaxy-by-galaxy ba- 
sis, albeit with som e assumptions about the reddening law. 
iBaldrv et al.l l|2008l ) show that these masses, determined by 
comparing star formation history templates based on the 



BC03 models to absorption features in the galaxy spectra, 
yield similar mass functions to other mass derivations using 
photometry, different population synthesis models, or differ- 
ent spectral features. 

In Figure [13] we show stellar mass functions (the num- 
ber of galaxies per Mpc^ per logarithmic stellar mass bin), 
separated into red and blue galaxies, for merger quenching, 
halo mass quenching, and SDSS galaxies. When dividing our 
simulated galaxy sample into red and blue, we do not use 
our dust prescription, so that Figure[T3]shows the mass func- 
tions for intrinsically blue and red galaxies. Given the small 
level of contamination of actively star-forming galaxies on 
the red sequence, this choice does not significantly impact 
the results. 

In our models, the red galaxy mass functions behave 
similarly to the luminosity functions, as expected since for 
these galaxies r-band absolute magnitude is a good tracer 
of stellar mass. The red mass functions do not have a sharp 
knee, and the slope for halo mass quenching is slightly 
steeper than that for merger quenching. In fact, above M* 
(the characteristic stellar mass in a Schechter function fit, 
« 10^° '' M©), the red galaxy mass functions in the two mod- 
els are remarkably similar. 

The blue galaxy mass functions reveal a crucial differ- 
ence between our halo mass and merger quenching models. 
Halo mass quenching produces a precipitous drop in the blue 
mass function above ~ IO'^'^'^Mq, a critical stellar mass 
corresponding to the critical halo mass 1O^^M0. Mer gers. 



Star Formation Quenching 17 



-2 



-3 



E 
CP 

o 



-5 



-6 




Mergers 

Holo Moss: 10' 

A SDSS data 



9.5 10.0 10.5 11.0 
Log(stellar moss) 



1 1.5 



Figure 13. Stellar mass functions for merger quenching (solid 
lines), halo mass quenching (dashed lines), and SDSS data (sym- 
bols). We split the galaxy populations into blue and red, as in- 
dicated by the colour of the lines. For merger quenching, the 
shape of the blue galaxy mass function tracks that of red galax- 
ies, roughly in agreement with observations. Halo mass quenching, 
however, yields a sharp cutoff in the blue galaxy mass function. 



in contrast, yield a blue stellar mass function whose shape 
traces that of the red stellar mass function (but with slightly 
lower normalization) for massive galaxies. The merger sce- 
nario markedly better reproduces the observed blue galaxy 
mass function. At M* > 10^^ Mq there are fewer blue galax- 
ies, but not zero blue galaxies. 

The overall conclusion from the stellar mass function 
comparison is the same as obtained from the luminosity 
function comparison: the merger quenching model produces 
red and blue stellar mass functions that broadly agree with 
data, but the halo mass quenching produces a sharp trunca- 
tion in the large blue galaxy population that is in disagree- 
ment with data. We emphasize that this truncation persists 
even in our variant quenching models which smear out the 
clumping associated with a halo mass cutoff (the hot frac- 
tion variant). This generic result is difficult to avoid given 
a relatively tight correlation between halo mass and stel- 
lar mass seen in the simulations. Unless this correlation has 
much greater scatter in nature than in our models, a sim- 
ple halo mass quenching model will have difficulty matching 
these data. 

Of course, the strict halo mass cutoff we adopt here rep- 
resents an oversimplification. Even when we tie quenching to 
the physically more relevant fraction of hot gas, we truncate 
star formation instantaneously and forever. Real galaxies are 
unlikely to behave this way. After a stable hot gas halo has 
formed and quenched the fuel supply, star formation may 
continue with an existing cold gas reservoir. Furthermore, 
passive galaxies may undergo mergers with gas-rich satel- 
lites, momentarily reinvigorating star formation and adding 
young, blue stars. 

Introducing a gigayear delay for the onset of quench- 
ing after a halo reaches the mass threshold, we still find 
no blue galaxies above IO^^'^Mq (while allowing more red 
gala:xies to grow stellar masses above 10^^ Mq). Such a time 
delay may not be much more realistic than our base model. 



but simulated galaxies around 10 Mq typically have gas 
fractions well below 0.5, so that they cannot grow by the re- 
quired factors ~ 10 using just their existing reservoirs of cold 
gas. This suggests that the steep cutoff in the blue galaxy 
mass function is robust. For pure halo mass quenching to be 
viable, we require a mechanism that chokes off the fuel sup- 
ply of most of the massive galaxies in hot haloes, but not all. 
A more detailed treatment of massive, star-forming galax- 
ies awaits fully hydrodynamic simulations that incorporate 
quenching mechanisms. 



6 DISCUSSION 

6.1 Origins of the red sequence 

Observations now su ggest red and dead galaxies have existed 
since at least z j:^ 2 ( Kriek et al.ll2008l : IWiUiams et~ai]|2009l : 



iBrammer et aL I l2009f r and the formation time for their stel- 
lar populations is quite old. When does the quenching take 
place in our models? 

Figure [14] shows the lookback time at which initial 
quenching occurred for each of the quenched galaxies in 
our two preferred quenching models. For the brightest red 
galaxies, quenching occurred at z ^ 3 in both models. In 
halo mass quenching, essentially zero galaxies brighter than 
r = —22 were quenched after z = 1.5, whereas in merger 
quenching a small number of bright galaxies came to the red 
sequence at late times. In halo mass quenching, red galaxies 
can only achieve the highest masses via mergers after mov- 
ing to the red sequence because the effective critical stellar 
mass (~ IO^'^'^Mq) is much smaller than the biggest galax- 
ies. In merger quenching, on the other hand, blue galaxies 
can attain quite high stellar masses and then move onto the 
red sequence as an already-bright galaxy at late times. 

Along the y-axis we show vertical histograms in tquonch, 
indicating the quenching rate in terms of the number of 
newly quenched galaxies as a function of cosmic time. The 
merger quenching rate peaks at z ~ 3, whereas a halo mass 
threshold quenches at a roughly constant rate over time. 
Here we caution that the peak in the merger quenching 
rate is sensitive to our minimum mass for merger quenching, 
which is set by the simulation resolution. A higher resolution 
simulation would likely move the peak to higher redshift, 
since small galaxies could quench earlier. In any case, the 
peak will remain as a distinguishing feature between merger 
and halo mass quenching. 

For a deeper physical understanding of our model galax- 
ies, we turn to Figure [TSl where we plot r-band stellar mass- 
to-light ratios, tyoung (the age of the youngest star parti- 
cle), mean stellar age, and mean stellar metallicity vs. stel- 
lar mass for our galaxies in our two preferred quenching 
models. We estimate mass-to-light ratios in solar units by 
dividing stellar mass by Lr, which is given by \og{Lr / du) = 
(r — C)/(— 2.5). Here r is the galaxy's r-ba nd absolute mag- 
nitude, C = -48.6 for AB magnitudes (|Okel Il974l ). and 
dv = d\[c/\^) is the approximate frequency width of the 
SDSS r-band filter. We take A = 6250 and d\ = 1500 
angstroms for the r-band. In both quenching models, we 
find a mass-to-light ratio that increases slowly with stellar 
mass, and is fairly tight above 10^°'^ Mq where most galax- 
ies have been quenched for some time. Halo mass quenching 
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Figure 14. Lookback time of initial quenching vs. 2 = r-band absolute magnitude for merger quenching and halo mass quenching. For 
mergers, tquench is the lookback time of the first major merger. For halo mass quenching, tquonch is the lookback time when the galaxy's 
dark matter halo first exceeded the critical halo mass of IQI^Mq. Along the t/-axis, we show histograms of tquench- The brightest galaxies 



'-'0 

all entered the red sequence at early times {z > 2). Merger quenching 
quenches galaxies with a rate nearly constant in time. 



quench • 

yields a peak quenching era at 2 3, whereas a halo mass limit 



displays a clump of M/L just below the quenching (stellar) 
mass, since these galaxies have recently undergone quench- 
ing. 

In the second row, we highlight a distinction in the age 
of the youngest stars between the merger quenching and halo 
mass quenching models. Massive galaxies in the halo mass 
quenching model lack any population of young stars: all the 
galaxies they accrete have old stellar populations, since even 
satellites that will be accreted later are quenched in massive 
halos at an early time. Galaxies quenched via mergers, how- 
ever, include trace populations of young stars. Even though 
mergers quenched the star formation in these massive galax- 
ies at z > 2, they obtain young stars via accretion of younger 
satellite galaxies. 

In the bottom panels of the figure, we examine ages and 
metallicities to understand why our models do not match 
the slope and normalization of the real red sequence. We 
include tracks showing the mean ages and stellar metallici- 
ties vs. stellar mass required to reproduce the observed red 
sequence. Specifically, on the plot of mean stellar age vs. 
stellar mass, we assume a uniform metallicity for all red se- 
quence galaxies, and ask: what mean ages would those galax- 
ies need to have to get the correct red sequence colours using 
the BC03 models? To answer this question, we create a grid 
of artificial single stellar populations (SSPs) with a variety 
of masses and randomly chosen metallicites and ages. This 
grid densely samples the region of the CMD where real red 
sequence galaxies lie. We then take our fit to the SDSS red 
sequence from Figured and for each absolute magnitude bin 
we identify all the artificial SSPs which fall within the scat- 
ter of the median m g — r and within 0.003 of the assumed 
metallicity. We compute the mean stellar mass and mean 
stellar age of these SSPs, and plot as a connected line in the 
figure. The plotted values have a scatter of ~1 Gyr. We fol- 
low an analogous procedure for the tracks in the metallicity 
vs. stellar mass plots. 

These tracks tell us what slope and normalization in 
age (metallicity) would be required to get the correct red 



sequence slope, assuming the red sequence has a constant 
metallicity (age). Of course, in real galaxies both age and 
metallicity may vary with stellar mass, but this is intended 
to illustrate the general trends. If age gradients along the red 
sequence are small, then we require a factor ~2 difference 
in metallicity from the brightest galaxies to those an order 
of magnitude less massive in order to reproduce the correct 
slope. If metallicities are roughly constant across the red 
sequence, then we require variation of 3-4 Gyrs between the 
bright and faint galaxies for typical galaxy ages. Finally, if a 
typical age is 10 Gyr, then the simulation metallicities would 
need to be increased by ~ x 2 in order to match the observed 
red sequence amplitude. We will discuss the implications of 
these trends in i]6.4l 

6.2 The blueness problem 

In all cases our models yield red sequences too blue by about 
0.1 magnitudes m g ~ r, an error of ~10% in the luminos- 
ity ratio between those bands. For these galaxies, colours 
are determined by matching each constituent star parti- 
cle's metallicity and age to grids from BC03 stellar popu- 
lation synthesis models. While the discrepancy might origi- 
nate with_jHi£ertontiesii^^ models 
(e.g. ICharlot et al.lll996l : IConrov et al.ll2009l '). we attempted 
to mitigate this poss ibility by using the same models as 
iBlanton etlii] l|2005bl ') for the VAGC. We also tried the up- 
dated Chariot & Bruzual models (graciously provided to 
us by S. Chariot) that employ an updated treatement of 
thermally pulsating (TP) AGB stars, but these showed very 
minor difference compared to BC03 for our red sequence 
galaxies since TP- AGB stars do not contribute significantly 
to the luminosity in SDSS bands. Hence we explored other 
possible reasons for the blueness problem. 

There are two ways to make a galaxy too blue: give it a 
younger age, or a lower metallicity. Perhaps our simulations 
simply have distributions of galaxy ages and metallicies that 
are wrong. As shown in Figure [TSl our brightest red galaxies 
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Figure 15. Redshift zero r-band mass-to-light ratios, age of the youngest star particle, mean stellar age, and mean stellar metallicity 
vs. stellar mass for galaxies in the merger quenching (left) and halo mass quenching (right) models. All panels include only red sequence 
galaxies. In the mean stellar age panels, lines represent the mean stellar age required to get the correct colours of the SDSS red sequence 
assuming a fixed metallicity along the red sequence of 1.2, 1.6, or 2.0 Zq (as labeled on the plot; see text), where Zq ^ 0.012. The 
metallicity plot shows analogous lines assuming a constant age along the red sequence of 6, 9, or 12 Gyrs. 



have mean stellar ages > 10 Gyrs. This leaves little room to 
maneuver, as galaxies cannot get much older given the age 
of the Universe. 

The story is different for metallicity. The tracks in Fig- 
ure [15] suggest that the metallicities of our galaxies are sys- 
tematically low. Very few galaxies show metallicities above 
0.02 (metal mass fraction), but BC03 models require that 
the most massive galaxies have metallicities significantly 
above this to explain their red colours. Why might these 
galaxies be unrealistically metal poor? 

As described in iFinlator fc Davj (|2008l ). metallicities 
in un-quenched galaxies are determined by an equilibrium 
between less enriched gas accreted from the IGM and gas 
enriched by ongoing star formation. Metallicities tend to rise 
as a galaxy grows in mass, and even at fixed mass galaxies at 



lower redshift are more metal rich (jPave et al.ll2007l ). Once 
quenching occurs in our models, the metallicities are frozen 
into the existing star particles. If our simulations quenched 
galaxies too early, the galaxies would not have built up their 
metallicities sufficiently. 

Another uncertainty is the adopted supernova metal 
yield s. We use a set of y i elds fr om various authors, most no- 
tablv lChieffi fc Limongil (|2004 ) yields for Type II supernovae 
that dominate the stellar metal budget of these galaxies. The 
yields are not all that well constrained, and it is possible that 
they are higher tha n we have assumed. We note that our 
simulations assume lAsplund et al] l|2005h abundances, for 
which solar metallicity is a metal mass fraction of 0.0126, as 
opposed to the 0.02 fraction that is assumed for solar metal- 
licity in the BC03 models. We assume that the metallicity 
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scale in the BC03 models represents an absolute scale (S. 
Chariot & A. Bressan, private communication), so that e.g. 
a simulated galaxy that has "solar" metallicity of 0.0126 is 
computed using the O.GZq BC03 templates. 

Our enrichment and galactic wind models further com- 
plicate the interpretation of the metal content of simulated 
stellar populations. In our simulations, star formation in- 
duces both metal enrichment and galactic winds. The winds 
suppress enrichment of subsequently-forming stars by re- 
moving metals from the ISM; this process is poorly con- 
trained observationally, and our model could be overly effi- 
cient at metal removal. We note, however, that the ejected 
metals in our outflow model match obs e rvatio ns of IGM en- 
richment (|Oppenheimer fc Pavel [ioog . |2009| ). so it is not 
clear that one can lower the ejected metals substantially. 
Since, however, most metals end up in the IGM particularly 
at early times (|Oppenh eimcr et al. 2009a), it would not take 
large changes to the outflow model in order to double the 
metallicity retained in galaxies. 

Our simulations do reproduce the stellar mass- 
:as phase m etallicity relation for galaxies at z ~ 2 
Finlator fc Dave 2008 ), suggesting that metallicities for 
star-forming galaxies in our simulation are generally cor- 
rect. However, systemat ic uncertainties a t the .2-0.3 dex 
level could impact t he iFinlator fc Davj l|2008l ) compari- 
son with data from lErb et al.l (|2006l ). These uncertain- 
ties emerge from th e calibration of metall icities derived 
from observation s (cf. [Kewlev fc Ellisonll2008l ). and because 
IFinlator fc Pavel M compared aU si=ed galaxies to 
UV-selected observations, which may have h igher specific 
star formation rates and lower metallicities (|Ellison et al.l 
I2OO8I ). Hence there is plausibly room for ~ x2 variations in 
the metallicities of galaxies. 

Finally, the under-enrichment could be related to the 
star formation history in the simulations. In particular, 
at least for the merger quenching, one expects that a 
starburst may be associ ated with the merger event (e.g. 
iMihos fc Hernauist|[l996h . We do not explicitly include such 
a starburst when we apply our quenching, and it could con- 
sume a significant amount of gas to create an enhanced 
metallicity. However, we note that we only quench galaxies 
typically well after the merger event, since we quench only 
at fixed time intervals of the simulation snapshots. While 
our simulations lack the resolution to track the detailed star 
formation history during the merger, the integral of the SFH 
is typically close to what would be obtained at much higher 
resolution, since it is primarily limited by the overall gas sup- 
ply. Hence while the resolution limits may modestly impact 
star-formation histories, it is unlikely that higher resolution 
can boost the metal production (or stellar production) by a 
factor of two. Indeed, a simple comparison of our base sim- 
ulation with another simulation with better mass resolution 
by a factor > 3 shows negligible differences in the mass- 
met allicity relation fo r star-forming galaxies at z = 1 (see 
also lDav6 et al.l l2007f ) . For resolved galaxies at fixed stel- 
lar mass, fractional differences in metal content between the 
two simulations are no larger than about 5 per cent, con- 
siderably smaller than the mean fractional scatter within a 
single simulation of 13 per cent. 

It is possible that there are some effects in the ob- 
servations that make the observed galaxies too red. De- 
riving reliable absolute magnitudes from apparent mag- 



nitudes requires several steps and important assump- 
tions. One must use some galaxy template t o calculate 
a galaxy's rodshift /distance and k-correction (iHogg et al.l 
,2002; Bla nton fc Rowei s 2002) from photometric data. In 
the VAGC, this is done by creating a suite of star forma- 
tion histories with corresponding spectra from BC03 mod- 
els, running a principle components analysis with the data 
to determine which star formation histories are most repre- 
sentative, then fitting each galaxy's photometric data to the 
principle components as templates. The complexity of this 
procedure makes it difficult to independently determine the 
reliability of the fits or w hether the chosen temp lates are 
sufficiently representative. iBlanton et al.l l|2005bf ) do esti- 
mate uncertainties in reported absolute magnitudes for each 
galaxy in the VAGC; for g and r, these are typically ~ 0.04 
magnitudes. Only systematic uncertainties much larger than 
these could give rise to the blueness problem, which seem 
unlikely. 

In summary, the blueness problem is likely to reflect 
some issue with the chemical enrichment level of simulated 
galaxies. The most obvious explanation is that the yields 
are somewhat too low, or that our outflow model is too ef- 
ficient at metal removal. In either case, while the overall 
metallicities are too low, the relative metallicities are basi- 
cally correct, which means that predictions for the slope and 
scatter of the red sequence are robust. However, it cannot 
be ruled out that this is reflective of a more subtle issue with 
the star formation or chemical enrichment histories of our 
simulated galaxies. 



6.3 Star formation histories 

By shutting off star formation in the largest galaxies, 
quenching mechanisms strikingly impact global star forma- 
tion. In our simulations without quenching, the most intense 
star formation occurs in the most massive galaxies, account- 
ing for significant fractions of the of global star-formation 
rate density. 

We show the effects of merger quenching and halo mass 
quenching on the star formation histories of galaxies of dif- 
ferent masses in Figure 1161 We first select galaxies from 
three stellar mass bins, 10^°, 10^ \ and 10^^ Mq, at z = in 
the simulation with no quenching model applied. The bins 
have a width of 0.1 dex in stellar mass. We then flnd the 
time each star in these galaxies was formed, translating to 
star-formation rate (in Mq/jt) by averaging over intervals 
of ~100 Myrs. Dividing by the number of galaxies in each 
mass bin, we get the average SFR per galaxy as a function 
of cosmic time. This methodology includes in the SFR any 
star-formation which occurred in small galaxies that later 
merged to form a massive galaxy that falls within one of our 
mass bins. 

For the quenching models, we track the same galaxies 
selected above from the no-quenching simulation, but we 
ignore the contribution of quenched stars particles to the 
SFR. We thus show the effects of quenching on speciflc col- 
lections of galaxies, rather than galaxies in given mass bins 
selected for each model (the quenching models do not pro- 
duce IO^'^Mq galaxies). 

Both merger quenching and halo mass quenching 
strongly suppress star formation in massive galaxies, while 
barely changing the average SFR in the lowest mass bin. 
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Figure 16. Average star formation rates (in Mq/jv) vs. cos- 
mic time for galaxies in 3 mass bins of width 0.1 dex: 10^", 10^^, 
and 10^^ Mq. We sort galaxies into bins by their 2 = stellar 
mass in the no-quenching simulation, then track the mass growth 
of those same galaxies through cosmic time in the no-quenching 
(black), merger quenching (green), and halo mass quenching (yel- 
low) cases. By construction, star formation within small galaxies 
that later merge to form a large galaxy is counted toward the av- 
erage SFR. Both quenching models strongly suppress star forma- 
tion in the two massive bins, but only slightly change the low-mass 
bin. Qualitatively, this behavior mimics the observed "downsiz- 
ing" of star formation. Halo mass quenching suppresses growth of 
the most massive galaxies more strongly than merger quenching 
because in halo mass quenching a massive galaxy's satellites will 
occupy the same halo, and thus be quenched, well before they 
merge. 



This effect natural l y leads to "downsizinR" with cosmic tim e 
ICowie et al.lll996l : [Heavens et alll2004l : ljuneau et aT]|2005h , 

where more star formation occurs in lower-mass galaxies at 
later times. Halo mass quenching suppresses SFRs more ef- 
ficiently than mergers because it quenches satellite galax- 
ies which later merge with the central galaxy. In contrast, 
merger quenching allows satellites to form stars until the fi- 
nal merger. Hence in principle, deriving accurate mean star 
formation histories from stacked galaxy samples in various 
mass bins can provide strong constraints on quenching mech- 
anisms. 



6.4 Successes and failures of quenching models 

The most basic result of this paper is that both merger 
quenching and halo mass quenching successfully recover 
qualitative aspects of observed local CMDs and LFs. These 
quenching models greatly improve upon the results of hydro- 
dynamic simulations that do not have any explicit quench- 
ing, which grossly overproduce the abundance of massive 
blue galaxies and fail to produce hardly any massive red 
galaxies. At absolute magnitudes —22.5 < r < —20, lumi- 
nosity function results from both models using our favored 
parameters fall within a factor of ~ 2 of observed luminos- 
ity functions from SDSS. The quenching models also yields 
a qualitatively reasonable colour-magnitude diagram when 
we apply a dust correction for star-forming galaxies. Qual- 
itatively, the red sequence grows as expected from observa- 



tions, with the most massive present-day galaxies truncating 
star formation at 2 > 2, and galaxies around L* continu- 
ing to evolve onto the red sequence to z = 0. This leads 
to mean stellar ages of massive red seq uence galaxies that 
are ro ughly consistent with observations (iGraves et al ]|2007l . 
I2OO9I ). These are notable successes that show that quenching 
associated with one or both of these mechanisms is on the 
right track towards understanding the evolution of massive 
galaxies. 

Models based on quenching hot mode accretion fail to 
produce a red sequence. This is mainly because even massive 
galaxies obtain a non-trivial amount of gas via cold mode 
accretion at the present day l|Keres et al.ll2009bl ). which is 
enough to keep massive galaxies too blue. Hence one must 
prevent the vast majority of cold gas from entering into mas- 
sive galaxies, be it pristine of recycled in a wind. Preventing 
hot and wind mode reaccretion produces a red sequence, 
but its LF is a power law with no obvious turnover at high 
masses. 

Despite their broad successes, the merger and halo mass 
quenching model still fail in subtle ways. Investigating these 
failures critically can help us to identify the underlying phys- 
ical processes that may be absent from our simulations. 

The main failures of the merger quenching model are (a) 
an excess of very bright red galaxies, (b) a shallow slope in 
the red sequence color-magnitude relation, and (c) a shallow 
slope and a less pronounced knee in the luminosity function 
(as well as the stellar mass function). 

The first problem results at least in part from post- 
quenching mergers. Quenched galaxies at the centres of mas- 
sive halos build up stellar mass by accreting satellites. We 
illustrate the importance of this stellar mass accretion in 
Figure 1171 which shows the z — stellar mass divided 
by the stellar mass at the time of quenching as a func- 
tion of z = r-band absolute magnitude. Since we show 
only quenched galaxies, post- quenching mergers drive all of 
the mass growth. In both merger quenching and halo mass 
quenching, the biggest galaxies have grown by factors > 10 
via mergers since the time they moved onto the red sequence. 
However, the merger quenching is generally larger at a given 
luminosity, and extends to factors up to 100. 

Too much mass is gained through this process, lead- 
ing to excessively bright galaxies: galaxies brighter than 
r = -22.5 (M. > IO^Mq) at z = in our model have 
grown in mass by a median factor of ^^3 since z — 1. This 
is in consistent with observational constraints (|Cool et al.l 
|2008| ) . Late-time mergers also explain the presence of young 
stars. Recent observations suggest no young populations in 
the b rightest red sequence galaxies (|Sanchez-Blazquez et al.l 
l2009l ). although up to ~ 30% of red sequence ga laxies do 
have ~ 1-3% of their mass i n stars < 1 Gyr old (|Yi et al.l 
I2OO5I : ISchawinski et al.ll2007l : iKavirai et all 120071 ). It is un- 
clear whether the presence and amount of young stars in our 
simulated quenched galaxies is inconsistent with these data; 
we leave a more detailed investigation for the future. 

By building up mass with small galaxies, post- 
quenching mergers may also contribute to the shallow slope 
of the red sequence colour-magnitude relation produced 
in the merger model. Stellar metallicities drive the slope 
of the red sequence, and the mass-metalli city relations 
(ITremonti et al. 1 12004 iFinlator fc Pavel |2008| ) tell us that 
massive star-forming galaxies are more metal-rich than their 
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Figure 17. Mass growth of quenched galaxies via mergers. The y-sxia is the log of the ratio of the stellar mass of quenched galaxies at 
z = to the stellar mass at the time of quenching. The upper envelope in merger quenching arises due to our resolution limit defined by 
a cut in stellar mass (see i|3.1l l. In both merger quenching and halo mass quenching, massive galaxies have obtained most of their stellar 
mass via mergers since they were quenched. 



less massive counterparts. If the mass-metallicity efltect dom- 
inates over the time evolution of metallicity at fixed mass, 
then a red and dead galaxy built up through accretion of 
many smaller galaxies will have a lower metallicity than that 
of a galaxy which grew just as massive through star forma- 
tion then moved to the red sequence. Any accreted satellites 
will have a lower metallicity, diluting the total metallicity of 
the galaxy. 

The shape of the luminosity function is more difficult 
to interpret. Its value in any luminosity bin is determined in 
a complicated way by the number of galaxies moving to the 
red sequence from the blue cloud, and the number of galax- 
ies moving along the red sequence. In principle, decreasing 
the mass accreted through satellites should steepen the red 
sequence slope, pushing the too-luminous galaxies to fainter 
r magnitudes, as well as keeping their metallicities higher. 

How might one lower the amount of satellite accretion? 
One possibility is that the issue is numerical. Simulations 
with low resolution tend to overmerge dynamical systems 
owing to stronger tidal forces on satellites whose spatial ex- 
tents are artificially enlarged by force softening. Another 
possibility is that much of the satellite mass i s actually go- 
ing in to diffuse halo (or intracluster) stars. IConrov et al.l 
l|2007h suggests that even fairly sizeable galaxies must be 
disrupted upon infall into clusters on a fairly short (~ 2 Gyr) 
timescale. Both trends suggest that our simulations proba- 
bly overestimate the growth of the central galaxy, at the ex- 
pense of intermediate-mass galaxies. Unfortunately, running 
higher resolution hydrodynamic simulations with sufficient 
volume to study very massive galaxies is not feasible, unless 
one performs constrained realizations of dense environments 
(in which case statistical comparisons become more compli- 
cated). 

Turning to the halo mass quenching scenario, the main 
failures are a) an excess of very bright red galaxies, b) es- 
sentially zero slope in the red sequence, c) no knee in the 
luminosity function, and d) a steep mass cut-off for blue 
galaxies. The first three problems are similar to those for 
merger quenching. Very bright red galaxies build up to an 



excess mass even though the satellites they accrete have been 
quenched. The problem is slightly less severe as in the merger 
quenching case, likely because satellite galaxies halt their 
star formation earlier. However, this trend also causes the 
slightly shallower red sequence slope as compared to merger 
quenching, because the accreted satellites are less enriched. 
When galaxies first move to the red sequence, their metal- 
licities are governed by the mass-metallicity relation. Since 
all red sequence galaxies arrive with roughly the same mass, 
they all start with roughly the same metallicity. Subsequent 
mergers between pairs of red sequence galaxies will then re- 
sult in a remnant of the same metallicity. This results in 
almost no trend with metallicity along the bright red se- 
quence (r < —22) in this model, as seen in Figure [T] which 
is inconsistent with the rather pronounced slope in the ob- 
served red sequence. 

The last problem appears more severe and intrinsic to 
the halo mass quenching scenario. The local universe clearly 
contains massive, bright blue galaxies; these are absent un- 
der this scenario. The difficulty is that the sharp halo mass 
cut produces a strong truncation of blue galaxies above a 
certain stellar mass, corresponding to the typical halo mass 
of truncation. Some blue galaxies are observed to have stellar 
masses well in excess of IQ^^Mq, which under this scenario 
would have to live in halos of mass < IO^^A/q. This is a 
remarkably high efficiency of star formation from the avail- 
able baryonic reservoir of such halos, well a bove what is in- 
ferred from observations at any halo mass (jMcGaueh et al.l 
l2009t ). This difficulty for halo mass-based quenching was 
also was found in a comparison of SAMs to SDSS data by 
iKimm et all ()2009( ). 

Overall, merger quenching appears to do a reasonable 
job of producing a red sequence as observed, and the dis- 
crepancies may be traceable to numerical issues. Of course 
a fully realistic model with quenching implemented dynam- 
ically within the simulation may alter the results, but the 
results here provide useful intuition and a starting point for 
more sophisticated models. 

How do our results compare to other galaxy evolu- 
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tion m odels ? The SAMs oflCroton et all (|2006l '). lBower et all 
(|2006l '). and ICattaneo et al I l|2006l ). which use broadly simi- 
lar quenching prescriptions to those we have adopted, appear 
better to match the observed color-magnitude relation of red 
sequence galaxies. These models have more freedom in the 
overall star formation histories of galaxies prior to quench- 
ing, since they contain free parameters to describe the gas 
dynamics. We also note that these models tend to populate 
a substantial red sequence without the additional quenchinq 
owing to hot haloes and/or AGN feedback, while our mod- 
els do not. A major difference in the gas content of haloes 
may result from our self-consistent treatment of outflowing 

? :as, which results in substantial re-accretion at late times 
Oppenheimer et al.l [ioOQh) as compared to tre atments of 
outflows in these SAMs. iDe Lucia et al.l (|201(]| ) point out 
how senstitive feedback prescriptions are to the gas content 
within haloes, so this could impact the nature and strength 
of quenching feedback. We are currently engaged in more de- 
tailed comparisons to SAMs, so we defer further discussion 
to future work. 



6.5 Generic problems in massive galaxy evolution 

There are several fairly generic difficulties in understanding 
quenching in the context of the observed galaxy popula- 
tion and its evolution. These are similar among all models 
of quenching we have examined, and hint at an underlying 
failure in our understanding or our current modeling tech- 
niques. 

One issue is that there is an intrinsic tension between 
producing a strong knee in the luminosity or stellar mass 
function, and having blue galaxies exist in substantial num- 
bers up to the highest stellar masses. The former would im- 
ply a very sharp quenching of star formation at a particular 
mass (stellar or halo), while the latter implies that quench- 
ing is a much more gradual function of mass. It is unclear 
how one can resolve this tension. Neither of our scenarios 
are able to do so, with merger quenching yielding large blue 
galaxies, and halo mass quenching producing a more pro- 
nounced truncation of the most massive galaxies. 

A second issue is that the most massive galaxies are ob- 
served to have gro wn by a s urprisin gly small amount o ver 
the last ~ 10 Gyr (|Cool et a l. 2008; Banerii et al]|2009l '). A 
consequence of this is that we overpredict the bright end 
of the luminosity function in both our quenching models; 
this is even the case when the satellites quench immediately 
upon entering a quenched halo (in disagreement with data). 
One way to avoid this would be to have fewer or less massive 
satellites, or to avoid having them fall into the central galaxy. 
Perhaps more detailed simulations will reveal that our cur- 
rent results are tainted by numerical issues such as a lack of 
resolution, but even SAMs have shown similar issues. Un- 
derstanding the dynamical evolution of a coUisionless stellar 
and dark matter halo seems to be an unresolved problem in 
the context of massive galaxy evolution since z ^ 1. 

A final issue is one of chemical enrichment, which is 
perhaps more subtle owing to current uncertainties in super- 
nova yields and outflows. Nevertheless, it appears that our 
current simulations have difficulty enriching massive galax- 
ies to observed levels, and perhaps more significantly, have 
trouble preferentially enriching massive red galaxies to levels 
approximately twice that of ~ L* systems. This has tradi- 



tionally been a difficulty of hierarchical models that assem- 
ble the massive end of the red sequence by dry mergers with 
lower-mass systems, as this will tend to decrease the aver- 
age stellar metallicities of the most massive system. Hence 
either massive galaxies must be preferentially enriched to 
begin with, which is seemingly contradictory to the ideas 
that they quench early on when galaxies of a given stellar 
mass have lower metallicity, or else the satellites they acquire 
must be quite metal-rich, which causes tension with trying 
to limit the amount of late-time growth. There is the further 
problem of understanding why a-enhanceme nt grows with 
galaxy velocity dispersion (|Graves et al.|[2007l ), which again 
seems to place some stringent conditions on massive g alaxy 
evolution (|Calura fc Mencill2009l : lArrigoni et"aLll2010n . Our 
current models do not reproduce this relationship, suggest- 
ing broader difficulties with the chemical enrichment histo- 
ries of galaxies in these models. Given the sensitivity of the 
exact shape and amplitude of the red sequence to metallicity, 
it is clear that a full understanding of massive galaxies will 
require understanding both their star formation and chem- 
ical enrichment histories. 



7 SUMMARY AND CONCLUSION 

By combining cosmological SPH simulations with simple 
post-processing prescriptions for the quenching of star for- 
mation, we test several proposed quenching mechanisms 
associated with major mergers, halo mass quenching, and 
quenching of hot and recycled wind accretion modes. We 
compared the galaxy populations resulting from these sce- 
narios with observations of the 2^0 universe from SDSS. 
With reasonable parameter choices, our merger and halo 
mass quenching models vastly improve upon existing simu- 
lations by suppressing the excessive growth of massive blue 
galaxies and creating a well-defined bimodality in galaxy col- 
ors. However, even these models fail to match observations 
in detail, providing clues towards key issues in the modeling 
of massive galaxy eveolution. 
Our main results are: 

• Merger quenching and halo mass quenching successfully 
produce a red sequence distinct from the blue cloud, with 
luminosity functions roughly consistent with those observed. 

• Quenching of hot mode accretion alone does not pro- 
duce a viable red sequence, showing that late time accre- 
tion continues via cold mode in massive galaxies. Further 
quenching wind mode accretion produces a red sequence, 
but a highly discepant luminosity function. 

• Our preferred merger quenching model rules out the re- 
formation of disks after gas rich mergers, because such re- 
formation leads to the build up of too many massive galaxies. 

• The halo mass quenching model is quite sensitive to the 
threshold halo mass: Mc ~ IO^^Mq yields the best match 
to the red galaxy luminosity function. 

• Both of these models yield an excess of bright red galax- 
ies due to mergers after the quenching process, and both 
yield somewhat too few galaxies around ~ L*; these prob- 
lems are more severe in the merger-based model. 

• Both models yield red sequences with too shallow slopes 
in the colour-magnitude diagram, likely due to near-constant 
mean stellar metallicity along the red sequence. 
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• In both models, the red sequence is too blue by ~ 0.1 
magnitudes, likely owing to simulated galaxies being too 
metal-poor by ~ x2 compared to real galaxies. 

• In both models, the brightest red sequence galaxies are 
quenched at redshifts > 2, in general agreement with ob- 
served estimates. 

We also identified some features that distinguish the 
merger quenching model from halo mass quenching: 

• Massive galaxies quenched via mergers include trace 
populations of young stars, whereas halo mass quenched 
galaxies do not. This is true only if satellites in addition 
to central galaxies are quenched in the halo mass scenario, 
but if this is not the case then the massive end of the LF is 
grossly overpopulated. 

• Merger quenching yields a population of very massive 
blue galaxies, with a mass function shape similar to that of 
red galaxies. Halo mass quenching creates a cutofT for blue 
galaxies at a stellar mass ~ IO^^'^Mq associated with the 
Mc = 10^^ Mq. 

• For resolved galaxies (M* > 3.4 x IO^Mq), the rate of 
galaxies quenched via mergers peaks at 2 ~ 3, whereas halo 
mass quenching occurs at an approximately constant rate. 

• The luminosity function for red sequence galaxies cre- 
ated via merger quenching has a shallower (bright end) slope 
than that for halo mass quenching. 

Overall, the merger quenching model seems to fare 
somewhat better than halo mass quenching, particularly in 
terms of the blue galaxy mass function, though it still has 
significant discrepancies versus data. The generic failures in 
both models in terms of the shallow red sequence slope, the 
lack of a pronounced knee at L* , and the excess of mas- 
sive galaxies reflect that there may be generic failures of 
the manner in which we quench galaxies here. It is possi- 
ble that some of the issues are numerical, but it could be 
that the quenching feedback processes actually impact sub- 
sequent galaxy evolution in a manner not accounted for by 
the post-processing technique we employ here. 

These results give us general intuition about how the 
red sequence is formed, as well as highlight the key issues 
that must be tackled by models. A concerted effort along 
both observational and theoretical fronts will be required 
to fully decipher the implications of these various trends. 
In future work, we plan to implement physical models for 
quenching within simulations dynamically. By identifying 
halos and merger events on the fly in our simulations, and 
quenching as appropriate, we can realistically track the feed- 
back from quenching that influences subsequent star forma- 
tion in galaxies. Given the evidence of the present study we 
favor quenching mechanisms involving major mergers, but 
this is unlikely to be the entire story. This work is but an 
early step towards understanding the origin and evolution of 
the fundamental bimodality in today's galaxy population. 
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